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[made. Two of four nights studied are nights with F s . 

The Perkins model is derived in a frame of reference moving with the 
velocity of the neutral wind; the model is transformed to the rest frame to 

facilitate comparison with data. Several data handling techniques are 
introduced. In particular, an integration interval that remains constant in 
length, but follows the vertical motion of the peak of the F layer is used 
to obtain the field integrated quantities of the Perkins model. 

It' was found that the Perkins model describes the general time behavior 
of the field integrated Pedersen conductivity if recombination is included, 
and the boundary flux is neglected. The numerical agreement is usually 
within a factor of two. It was found that the flux data is inconsistent with 
the rest of the data. Errors on the order of 30 m/s in the individual 
velocity measurements are believed to be present to account for the dif- 
ferences . 

The Perkins model is able to predict reasonable growth rates of the 
development of F g ; however, a direct comparison between the growth rates ■ 
predicted by the model and actual onset times is not possible because, of the 
lack of data on the eastward neutral wind. 

The linear growth rate of the Perkins model is extended by including' 
recombination, and a short wavelength damping term and coupling to the E 
.region. It is found that, for wavelength perturbations several kilometers 
or longer, the E region can short circuit instabilities in the F region and 
reduce the growth rate . 
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ABSTRACT 


The current models of mid- latitude F are studied* The 

s 

assumptions and derivations of the Reid model, the Scannapieco 
model, and the Perkins model are presented in detail* 

Incoherent- scatter data of the density profiles and velocity 
profiles were obtained from the Arecibo Observatory in order that 
the models could be evaluated on the basis of experimental data. 

Initial studies indicated that the Perkins model was most representa- 
tive of the data from Arecibo, so a detailed comparison of the 
predictions of the Perkins model and the data was made. Two of 

four nights studied are nights with F . 

s 

The Perkins model is derived in a frame of reference moving 
with the velocity of the neutral wind; the model is transformed 
to the rest frame to facilitate comparison with data. Several data 
handling techniques are introduced. In particular, an integration 
interval that remains constant in length, but follows the vertical 
motion of the peak of the F layer is used to obtain the field integrated 
quantities of the Perkins model. 

It was found that the Perkins model describes the general time 
behavior of the field integrated Pedersen conductivity if recombination 
is included, and the boundary flux is neglected. The numerical 
agreement is usually within a factor of two. It was found that the 
flux data is inconsistent with the rest of the data. Errors on the 
order of 30 m/s in the individual velocity measurements are believed 
to be present to account for the differences. 
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The Perkins model is able to predict reasonable growth rates 


of the development of F ; however, a direct comparison between the 

s 

growth rates predicted by the model and actual onset times is not 
possible because of the lack of data on the eastward neutral wind. 

The linear growth rate of the Perkins model is extended by 
including recombination, and a short wavelength damping term and 
coupling to the E region. It is found that, for wavelength perturbations 
several kilometers or longer, the E region can short circuit instabili- 
ties in the F region and reduce the growth rate. 
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CHAPTER I 


INTRODUCTION 

Spread (F g ) is caused by the presence of plasma density 
irregularities aligned along the earth's magnetic field (B) in the 
F region of the ionosphere. The term "spread F" is derived from 
the spreading of the returned signal on an ionogram. An ionogram, 
which is a measurement produced by reflecting radio waves from 
the ionosphere, normally would have two discrete traces of the 
returned signal from the F region indicating the reflection height 
vs, frequency. During F , the field aligned irregularities cause a 
smearing, or spreading of these traces. The smearing can either 
appear as "range spreading", when there are many reflection heights 
for the same frequency, or "frequency spreading", when waves of 
different frequencies are reflected at the same height. Two ionograms 
are shown in Figure 1; the ordinate is the virtual height of reflection 
of the signal, and the abscissa is the frequency of the signal. One 
shows normal traces from the undisturbed ionosphere, and the other 
shows the spreading of the signal representative of F g , 

The density irregularities that cause F g can be grouped into 
two general classes: very large patches that extend horizontally for 
several hundred kilometers, and smaller individually aligned 
irregularities with sizes on the order of kilometers or less. It is 
likely that the large patches are composed of many of the smaller 
ones. The irregularities are seen from the bottom of the F layer to 
around 1000 km (Herman, 1966). 
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The occurrence of F is negatively correlated with magnetic 

s 

activity at the equator and positively correlated at the poles. In 
general, F appears more during spring and fall than summer and 

s i 

.winter, and is rarely observed between 20° and 40° latitude {Herman, 

1966 ). • 

Recently, Zinchenko (1976) studied the probability of occurrence 

of F at Arecibo Observatory as a function of season, time of day, 

and solar activity. He has shown that the Arecibo Observatory, is 

atypical of some of the general characteristics described above. 

In particular, Zinchenko (1976) found that Arecibo, which is at 30° 

latitude and is not expected to see F , observes F g quite often, 

especially in the winter and .summer months. In fact, during 

December and January, the probability of observing F on any given 

s 

night approaches 90 per cent. 

There is a great deal of data on F g occurring at equatorial 
latitudes. (Farley, et al„ , 1970; McClure and Woodman, 1972; 

Balsley et al. , 1972; Farley, 1974). Farley et al. , (1970) list some 
of the important characteristics of the irregularities that produce 
equatorial F g : 1) the irregularities usually appear when the F layer 

is high, and moving up. They may persist for an hour or more as the 
layer moves down, and then are quenched; 2) the irregularities are 
most often seen on the bottomside of the F layer; 3) they usually 
appear a few hours after sunset, but may be delayed if the upward 
motion of the layer is delayed; and 4) they are never seen during the 
day, even during a total solar eclipse. 
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The data base for mid -latitude F is not as extensive as that for 
, s 

equatorial F g . The important works are those of Zinchenko (1976), 
describing probabilities of occurrence of F g at Arecibo, and Kelley 
(1'972) and Dyson et al. , (1974) that show that F g occurs simultaneously 
with electrostatic turbulence. 

A number of theoretical models have been proposed to account 
for equatorial and mid-latitude F g . Those of Hudson and Kennel (1975) 
and Williams and Weinstock (1970) are applicable only, to equatorial 
latitudes. Reid (1968) and Cunnold (1969) proposed an E x_B drift 
. instability that is applicable at'all latitudes, and Perkins (1973) 
and Scannapieco et al„, (1975) propose detailed models that are 
'applicable to mid- latitudes only. Each of these models make certain 
assumptions about ionospheric conditions. Each postulates a definite 
set of onset parameters that must be satisfied before F g will occur. 
However, a comparison of the models and their assumptions with the 

available data has not been made. This, thesis will partially fill that 

* 

void, by comparing’ the current models of mid-latitude F g with data 
from the Arecibo Observatory. 

The derivations of the Perkins model, the Scannapieco model, 
and the Reid and Cunnold models are presented in Chapter II. The 
data that are available from the Arecibo Observatory are- introduced 
in Chapter III, Several data handling techniques, as well as methods 
of obtaining, ionospheric parameters from the raw data are also 
presented .in Chapter III. Also, the Perkins model is extended, and 
the -equilibrium conditions predicted by the extended model are 
compared with, the data from Arecibo. 
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In Chapter IV Lt is shown that the data available from the 

Arecibo Observatory are insufficient to calculate growth rates predicted 

* ! 

by the Perkins and the Scannapieco models,, However, it is shown that 
the models do have the capability of predicting reasonable (i.e„, 
physical) growth rates. The Perkins model of F g .is extended by 
including recombination, E region coupling, and short wavelength 
damping due to diffusion. 

The summary and conclusions are presented in Chapter V. 



CHAPTER H 


MODELS 

In. this chapter, three mid- latitude models will be discussed. 
These include the Perkins (1973) model,, the model of Scannapieco 
et al. (1975), and the Reid (1968) model. The Cunnold ,(1 969) model 
is similar to the Reid ,(1968) model and will not be discussed 
separately. The three models chosen .are the .most current and 
complete models relevant to mid- latitude F g . 

The derivations of the models are discussed in.fairly complete 
detail to enable the reader to better understand them. They also set 
the stage for the extensions that will be added. in later chapters. 

2.1 The Perkins Model. 

The basic premise of the Perkins model is that the nighttime 

F region is supported against gravity, by electrodynamic forces 

(j_xJB) and the neutral wind (Perkins, 1973). The current j_ is the 

Pedersen current, which is parallel to- the electric field (E) and 

perpendicular to the magnetic field .(B). The Hall current, which, is 

perpendicular to JE and_B, is assumed negligible, in- the F region 

2 2 

because it is of order (v / w .), where v . is the ion-neutral collision 

m c i in 

frequency and w . is the ion cyclotron frequency. The ratio v . . 

is much less than unity in the F region. The Perkins model also 
assumes that: (1) the ions and electrons are- isothermal; (2) 
recombination and ion drag are negligible because the time scales for 
these processes at night is on-, the order of hours, which, is much 
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longer than the time scale for the development of F ; (3) the magnetic 

s 

field is uniform; (4) the neutral atmosphere consists of a single species 
distributed exponentially; in altitude with a constant scale height, H; 
and (5) E region contributions, to the integrated conductivities are 
neglected. 

Two of these assumptions, the neglect of recombination and 

the neglect of E region contributions,- will be examined in later 

chapters. The neglect of recombination is not valid when, the F layer 

is low; the time scales then approach 15-30 minutes, which, is of the 

same order as the time scale for the development of F g . Perkins 

(1973) neglected the E region contributions to the integrated 

conductivity because he assumed that, the much greater extent of 

the F region would overcome the greater conductivity per particle- in 

the E region when 'integrated quantities were calculated. Zinchenko 

(1976) found that this as sumption, is not generally valid. 

The equations of motion are written. in a drift frame moving 

with the neutral wind (V ), which is assumed constant. The coordinate 

— n 

system that will be used exclusively in this thesis is shown in 
Figure 2. The view is toward the east in the northern. hemisphere. 

A /v 

The vectors labelled h and 1 point towards the vertical and northward 

A A 

directions respectively, while those labelled z andx.a-re parallel 
and perpendicular to the magnetic field respectively. The magnetic 
field is pointing downward, with a dip angle D from the northward 
direction. The equations of motion are the following: 
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ions: 


v. x B 

-TV n + ne — — 

— c 


ne $ + n M g 


-Mn v. v. = 0 
—i in 


electrons: 

v x _B 

-T V_ n - ne - - - e - + neV <5 = 0 

The view is toward the east in the northern hemisphere. The 
A /\ 

vectors labelled h and' 1 point towards the vertical and northward 

A A 

directions respectively, while those labelled z.andxare parallel 
and perpendicular to the magnetic field respectively. 1 The magnetic 
field is pointing downward, with a dip angle D from the northward - 
direction. The effect of collisions on electro.n motion has been 
neglected reflecting the fact that the electrons are tied quite strongly 
to the field lines in the F region. Gravity has been neglected 
because of the low mass of electrons. 

By solving -the electron equation of motion for <5 the 
potential can be separated into two parts. 


$ = <p-(x, y) + — In -3- 
w e n 

o 

( . 

The quantity n^ is a reference density. The equations 'of motion - 
may now be written. 
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ions: 

v . x B 

-2 T _V n + ne ~r - ne V_ <p { 

+ n'M g - M n v. v. = 0 (2. 1) 

° — i m 

electrons: 

v x _B 

^ + _V 0 £ 0 (2.2) 

Because the Perkins model is developed in a frame of 
reference moving with the neutral wind, the electric field (V <$ as 
written inEquations 2.1 and 2.2 is actually a resultant electric field 
that implicitly contains the neutral wind. It is of the form shown 
below: 


-RES - + 




where _E is the electric field that would be measured from a stationary 
reference frame, and is the neutral wind velocity (assumed constant). 
This point is quite important when one desires to compare .the Perkins 
model with data; electric field measurements are usually made from 
the ground, and thus are in the stationary frame. The parallel and 
perpendicular ion velocities and perpendicular electron velocity are 
obtained from Equations 2. 1 and 2, 2. 

g sin D _ 2T . _I_ 3 n 


v . 
zi 


n 3'z 


(2.3) 
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V 


i-i 


c /\ 2T 

- g- V 0 x * - — 


c 

B 


„ A 

V nx z 


n 


+ 


- $k* z 

to . 
ci 


in 


10 


ci 


# V, 0 + ^ T 


c 


^o. 11 


n 


Su 


w . 

Cl 


2 2 

+ order (v. /to .) 

v in ci 


(2.4) 


— e 


c „ , A 

— _V 0 x z 


(2.5) 


The ion continuity equation is 


+ V n 
8 t — ^ 


A 

3* z . 


to . 

Cl 


B 


V 


+ 


9 | ng sin D 2 T 9 n 

M v. 9 z 


9 z 


v. 

in 


- 0 


in 


( 2 . 6 ) 


Terms of-first order in v.^/to^ have been neglected in Equation 2.6 • 
because they are much smaller than the zero order terms. However, 
in the calculation of the transverse electrical cur.rent, the. zero . 
order terms cancel out, so. the first order terms. must be retained. 
This procedure of taking Equations 2,6 and 2.7 to-.two different orders 
in v^/oo^ leads to a slight inconsistency, which will be pointed out in 
Chapter III. 

The electric current perpendicular to the magnetic field (j^) 
is obtained from Equations 2.4 and 2.5: 
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Jn = 


2Tc v 
B — 


A 

n x z 


+ ne 


JL* 


A 

z 


CO 


Cl 


ec 

__ 


v. 

in. 


co 


n Zx 0 


ci 


2Tc- 

B 


v. 

in 


V n + 


co 


ci 



v. 

in 

co • 


ci 


n .&- 


(2.7) 


Two moments of the ion density distribution are introduced: 

N(x,y) = I nd z 

I it. ndz 

^ ci 

It is assumed by Perkins (1973) that the limits of integration can be 
chosen so that the density and derivative of the density vanish at the 
limits., Thi's is equivalent to assuming that the flux vanishes at the 
endpoints* This assumption does not represent the true F layer, as. 
will be shown. in the next chapter. The quantity N is the field 
integrated ion density, or the plasma content of a field tube, and S 
is the field integrated F region- Pedersen conductivity. 

In the absence of ionization and recombination no net curr.ent 
will enter or leave a flux tube. Consequently, 


5 . x “ y L dz = 0 


( 2 . 8 ) 
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This is really a statement of plasma quas i- neutrality, as can be seen 
by subtracting the electron continuity equation from the ion continuity 
equation, and integrating. By substituting Equation 2.7 into 2.8, and 
using the moment definitions of S and N, one obtains: 


V^SV^. <p - 


0N 
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g e cos D + _2T_ v 2 
(x> . e *• 


8 2 


8 x 


cos D 


2T 

eH 
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(2.9) 


It has been assumed that the ion-neutral collision frequency has the 
same height distribution as the neutral atmosphere. 


v. 

in 


v. 


in 0 


-h/H 
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The ion continuity equation, Equation 2. 6, is integrated along 
a field line to obtain the equation describing the time evolution of 
the plasma content of a flux tube: 


8 N 
8 t 



cos D 


8 N 

9 y 
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-|r" ^ x % - 0 
B — 

( 2 . 10 ) 


To complete the description, one also needs an equation for the 

time evolution of the integrated Pedersen conductivity. This equation 

v. 

is obtained by multiplying Equation 2.6 by — ^ , and integrating 

^ci 

along z: 
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( 2 . 11 ’) 


Equations 2.9, 2. 10 and 2. 11 form, a complete set of "equations for 

the description of N, 2 , and <p° 

The equilibrium state, i.e., the state before F , is one in 

s 

which there are no gradients in'N and 2 transverse to B, and a 
transverse electric field is present. The equilibrium solution of 
Equation 2. 11 is a requirement that there be an eastward Pedersen 
current to support the F layer against gravity: 
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co . cos D 
ci 


{ 2 . 12 ) 


The F g model is based on the fact that this equilibrium is unstable to 

small perturbations. This is shown .by linearizing Equations 2, 9, 

2 . 10 , and 2.11 about the equilibrium, and assuming a solution of the 
ik ° x *f v t 

form e— for the perturbed quantities in a frame of reference 

moving with the velocity 


. g x B 


ci 


IT 1 *0 x * 


No perturbations are allowed in N because there is no way to change the 
plasma content of a flux tube in the absence of recombination or 


ionization. 
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The complex growth rate y obtained by the perturbation 
analysis is: 
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(2* 13) 


The quantity is the neutral mass, the angles u and 0 are shown in 
Figure 3a, and 
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It should be re-emphasized that the electric fields in Equations 

2. 12 and 2. 13 are resultant electric fields as shown in Figure 3b. A 

positive value of real y is required for growth; a negative value means 

the perturbation will be damped. 

The first term in brackets in Equation 2. 13 will dominate 

because v. << oo . . The maximum growth rate is obtained when 
in ci & 

a = 0/2. 



X (NORTH) 





Y (EAST) 

SOUTH 

3a: Orientation of the E and k vectors. 
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Figure 3b: The production of the resultant electric field from, 
the electric field and the neutralvwind. 
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Y MAX 
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H < v. - > 
in 
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sin D 


sin 0/2 


cos 0 


(2.14) 


Equation 2.14 predicts, that F‘ will be- more .likely to occur when "the 

s 

layer is high, because 1 / <v increases with height. This is. in 

agreement with observation (Herman, 1966). It also predicts a 

maximum growth when there is- a large north-south component of the 

equilibrium resultant electric field; cand the orientation. of the striations, 

which is the direction perpendicular tok, will be more north- south 

than east- west. It should be noted that there is no' wave length 

dependence in the Perkins growth rate. It will be showndn Chapter IV 

that the wavelength dependence, is contained in ther terms of order 
2 2 

v . /co . that were neglected in.the Perkins derivation, 
in • ci ° 

This instability can.be described physically with the help of 
Figure 4. The ionosphere has been perturbed. in such a way that 
there are rising and falling sheets of ionization aligned along the 
field lines. In. region 2 the F layer has fallen so that 2^ has increased 
over the equilibrium value 2^, in region. 1. Because V ° _J = 0 and 

V_ x E = 0, the normal component of the current, which is assumed' 
to be a Pedersen current, must be continuous, and the tangential 
component of the electric field must be continuous. Since J - 2E 
an d > Sj then 2 > and > E^. Thus, the supporting 
eastward current has been decreased in region 2 as shown in Figure 4, 
and the F layer will fall even further, causing 2 2 to become even 
larger. Because there is no damping .mechanism, the process-will 
continue and the system as shown is unstable. 
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2. 2 Scannapieco Model 

The Scannapieco model, as described in two recent papers by- 
McDonald et al., (1975), and Scannapieco et al., (1975), starts with 
the same flux tube equations as. the Perkins model (Equations 2, 9, 

2,10 and 2.11), However, by the assumption of a different equilibrium 
state, a totally different model of F g is obtained. It is of the E x B 
gradient instability type, as will be shown. 

What makes this model different from the Perkins model is the 
assumption of an east-west gradient in the field integrated Pedersen 
conductivity. The equilibrium is studied in a reference frame -moving 
eastward with the velocity: 
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The equilibrium condition on the Pedersen current is obtained from 
Equation 2.11: 
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(2. 15) 


This condition must also be consistent with Equations 2. 9- and 2. 10. 
The following is obtained by integrating Equation 2. 9-. in y: 
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(2. 16) 
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By comparing Equations 2* 15 and 2„ 16, one obtains; 


CONSTANT = 


" N o e g 
cos D co 


c i 


It is not inconsistent with Equation 2. 10 to have N q constant so there 
is internal consistency among the three equations. Actually, the 
assumption that N q is a constant can be justified by the data, at 
least for time scales on the order of an hour. This will be demonstrated 
in Chapter III. 

The growth rate analysis begins with the assumption -that all 
perturbed quantities have the form ex Y . Propagation of the 

irregularities is only. allowed in the x direction because of the gradient 
in S in the y direction. , A perturbation propagating in the y direction 
can not be allowed in this normal mode analysis unless the wavelength 
of the perturbation is much less than the characteristic length of 
the equilibrium gradient. - 

By linearizing Equations 2.9, 2, 10 and 2.11, and substituting 
the assumed solution for perturbed quantities, one obtains. the following 
expression for the -real part of y , which is the growth rates 
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The final form is obtained by substituting Equation 2. 15. for — ~ 
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,(2. 17) 


If there is no- Pedersen gradient, a pure damping solution, is' obtained,' 
which is the same as the Perkins model result for propagation, in the 
x direction (c.f. Equation 2 . 13). 

This model of F is really, of the E x_B_ gradient instability class, 

as can be demonstrated with the help of Figure 5. The ionosphere has 

been perturbed such that alternating sheets of Pedersen. conductivity 

exist, oriented purely- in the east- west direction. In region 2, the F 

layer has risen, so < Sj. The x component of the electric field 

will cause E x B drifts towards the east, as shown in Figure 5. 

Because the current normal to the striations must be continuous, 

S. E , -= E Since 2L < 2, , then E ~ >E .. . Thus, the eastward 
1 nl 2 n2 ^1 n2 nl 

drift velocity in region 2 will be greater than that in region 1. Region 
2 is a region of depleted Pedersen conductivity moving rapidly into 
larger ambient conductivity<;because of the equilibrium gradient. This 
makes the value of 6 S/S o even greater as time goes on, and hence it is 
unstable. Region 1 remains stable. 




Figure 5: Physical model of the Scannapieco instability 







The gravity term in’ Equation* 2. 17 would normally be expected 
to be small. However, it can contribute to the- ins lability because -the 
x _B drift would cause the whole layer to move eastward into the 
ambient gradient. 

The constant damping term in Equation. 2. 17 represents the 
linear increase of 2 as the layer falls down the field lines. As the 
Pedersen conductivity in region 2 becomes smaller due to the rising 
of the F layer in that region, the tangential electric field must be 
continuous. This causes the tangential current, which- is the 
supporting eastward current, to decrease. The layer falls and 2 
. increases-, thus contributing a damping mechanism. 

2.3 Reid Model 

The Reid model, presented in a paper by Reid (1968), and 
extended later by Cunnold (1969), is totally different from the two 
previous models discussed in that it is not a flux tube model. No 
quantities are integrated. along the field lines; it is thus a local 
model and its predictions are a function of altitude. 

The Reid model is of the _E x_B instability type, and depends on 
the fact that the transverse mobilities are different for the electrons 
and ions. .It can be described with the help of Figure 6. The solid 
horizontal line represents the unperturbed electron (and ion) density. 

If the density is now perturbed, as in the sinusoidal curve, the whole 
plasma will move to the left due to the _E_x B_ drift. The greater Hall 
mobility of the electrons will enable them to move further to the left 
than the ions. This is shown on the figure by using two sinusoidal 




Figure 6: Physical model of the Reid instability. 
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curves slightly out of phase. As shown, space changes will develop 

I 

causing perturbed electric fields (E ) which will in turn cause 
e' x_B_ drifts. The region of excess ionization will move into a region 
of less ambient ionization, as shown by the direction of the ambient 
density gradient, while the region of depleted ionization will move 
into a region of greater ambient density. The process will continue 
in the absence of a damping mechanism causing Sn/n Q to further 
increase. The same effect could occur with an electric field directed 
to the right since the Pedersen mobilities are different for electrons 
and ions. 

• The equations used in the Reid model are the equation of 
mo tion 


m v 

- V n T + n- E + j x.B - S Sn 1=0 

— s s — — s — q -*s 

^-s 


and the continuity equation: 
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9T“ 


-f V • n-.v 
— s— s 
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The subscript "s" represents the species, i, e, "e" for electrons and 
"i" for ions. Thus there are four equations. 

The equation of motion is solved for the particle flux which is 
then substituted into the continuity equation. The two resulting 
equations are linearized, and a solution of the form 

-x cos D -i(o)t - k y) 

H 6 ? 
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exp 



27 


is assumed for all perturbed, quantities. The quantity is the 
plasma scale height. „ The Reid model only allows propagation in the 
eastward direction. 

The growth rate is found by assuming that the wavelength of 
the perturbation is less than a few kilometers, and H^.is greater than 
5 km, i.e. kH > > 1. By substituting the assumed solution in the 
continuity equations for electrons and ions, one obtains, the growth 
rate: 
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Y * = H ^ 1 + ( i *)" 1 

i zt 

i = Ls the scale height of . 

The quantities D p> D^., jjl , and ^ are the Pedersen and Hall 
components of the diffusion coefficients and mobilities, respectively, 
and are defined as 


D = D 


. ±2 

(v ) 


o . ± 2 . =fc 2 

(v ) + (u c ) 


H 


D* 
' o 


± ± 
W V 
c 


, ±.2 t . ±2 

(v ) + (co c ) 


The mobilities have the same form, and 
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The requirement for the onset of the instability is that o>. >- 0. 

2 1 
If ky ^ then the onset condition is 


E ox (K 1 L 4 cosD + Kj K 3 cos D) - E q ^ (Kj L 2 cos D) 


(K 2 K z .k^ + Lj E 3 cos ^ D) 
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with a minimum wavelength given by: 


X MIN < 


E cos D - E 
oy ox 
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If the equilibrium electric field (E ) is increased, the minimum 
wavelength and growth time are decreased. The growth times are 
also dependent on latitude, becoming infinite at the poles where 
D = 90°„ 

Reid (1968) proposed that irregularities are produced by this 

instability in the E region, where the differences in the mobilities are 

the greatest, and are then transmitted to the F region via the field 

lines causing F g „ No analysis was performed to determine the 

conditions for such propagation The instability is damped at high 

latitudes because of the dip angle D, but there are increased electric 

fields at high latitudes to compensate,, At mid- latitudes, there are 

no large fields, and the instability is damped, thus explaining the 

absence of F between 20° and 40° latitude. The correlation with 
s 

magnetic activity is explained also. During magnetically active times 
at high latitudes there is much uneven particle precipitation, causing 
increased density gradients, and stronger electric fields. These two 
conditions will enhance the instability, thereby explaining the positive 
correlation. At the equator, there is a strong return current during 
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magnetically active times that will oppose the normal post- sunset 
current. The electric field is decreased, explaining the negative 
correlation. 

Cunnold (1 969)' extended the Reid model by including electron 
ion collisions, losses due to recombination, and propagation in an 
arbitrary direction. The positive part of the growth rate obtained 
by Cunnold (1969) is 


k K 


(V . - V ) 
— oi — oe 
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V = drift velocity due to the zero order electric field.. 

— o 

This term must overcome any damping effects (e.g. diffusion) to 
have a growing solution. 

Perkins (1973) makes the following criticisms of the Reid 
model: 1) the growth rate changes sign as a function of altitude; 

2) the question of whether there is a bona-fide equilibrium about 
which to perform stability calculations was never- brought up;, and 

3) the vertical dependence of the perturbation amplitude was treated 
in an ad-hoc manner. 


2.4 Statement of the Problem 

•Each of the three models discussed in the previous sections 

describe an instability that could produce F on an ionogram. 

s 
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However, none of the models has been adequately tested by comparison 
with the data that are currently available on the nighttime ionosphere. 
The Perkins model is chosen for detailed comparison with the data 
for the reasons described below. 

The Reid model is an E region model. The instability depends 

on the fact that the mobilities for motion perpendicular to the magnetic 

field are different for electrons and ions. Because the collision 

frequencies are so large in the E region, transport across field lines 

is easy, and the mechanism described by the Reid model could 

easily occur. In fact, Reid predicted growth times on the order of 

several minutes using typical ionospheric parameters for the E region. 

Ions and electrons in the F region, however, are much more tightly 

bound to the field lines because of their lower collision frequencies.' 

Thus, the _E^x _B_ instability proposed, by Rei'd would not be effective - 

in the F region. This has been checked by the data introduced in 

Chapter III and it was found that even on the bottomside of the F region, 

where the gradients and collision frequencies are largest, the growth 

times predicted by the Reid model are on the order of hours. This ' 

is not physically meaningful in the case of F , which can appear in a 

s 

matter of one half hour. 

Reid (1968) proposed that his model could produce irregularities 
in the E region, and these irregularities could then produce perturba- 
tions in the F region. However, recent work by Zinchenko and Nisbet 
(1976) shows that the perturbations of wavelengths shorter than 5-10 km 
will not be transmitted from the E to the F region. Consequently, the 
Reid model has been ruled out as a viable means of explaining F . 
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The Scannapieco model, which allows equilibrium gradients in 
the Pedersen conductivity is essentially an extension of the Perkins 
model. However, for the four nights studied in this work, the 
existence of the horizontal gradient is not found, except at sunset. 
Scannapieco et al. , (1975) state that the gradient at sunset would 
remain into the night because it is an equilibrium solution. The data 
from Arecibo Observatory which is discussed in Chapter III does not 
support this contention. In fact, the east-west horizontal gradients 
seen are quite small. The growth times predicted by the Scannapieco 
model using the data available are even longer than those predicted by 
the Reid model. However, it should be pointed out that it is difficult 
to acquire data on the horizontal structure of the Pedersen conductivity, 
and the data that was used to obtain the growth times mentioned above 
may not be accurate enough to make a definitive statement about the 
existence of a gradient. Discussion of this point will be continued in 
Chapter IV. 

The Perkins model was chosen for detailed study because it 

has the capability of predicting reasonable growth times (10-30 minutes) 

and the equilibrium state described by the model in which the nighttime 

F layer is supported by JE_x _B_forces and neutral winds is generally 

accepted as true (Behnke and Kohl, 1974; Harper, 1973; Behnke and 

Harper, 1973). The Perkins model is a flux tube model, reflecting 

both the field aligned nature of the striations associated with F and 

s 

the fact that plasma particles in the F region are tightly bound to the 


field lines. 
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A major gap in the literature on F is a comparison between 

s 

the models of F g and the available data. This thesis will attempt 

to fill this gap by comparing the Perkins model of F with data from 

s 

the Arecibo Observatory in- Puerto Rico. The Perkins model is also 
extended based on the findings of the data analysis. 



CHAPTER III 


EQUILIBRIUM ANALYSIS 

3. I Derivation of the Equations in the Laboratory Frame 

Certain changes must be made to, the Perkins model before it 

can be compared with the data. The Perkins model, as well as the 

other models discussed in Chapter II, was derived in a frame of 

reference that is moving with the constant velocity of the neutral wind. 

However, the data is obtained in the laboratory frame. The transforma- 

]B A * 

tion to this frame introduces an extra term, -vV.^ x z, that is added 
to the electric field. As was indicated in Chapter II, a resultant 
electric field is produced, and is of the form shown below: 


C TT 

B -RES 


- * 


+ V 
+ -N 


A 

X z 


The electric field terms that appear in the models are actually 
resultant electric fields of the form shown above. The authors of 
these models are not neglecting the effect of the neutral wind, but 
are implicitly including it in the electric field terms. It is also 
assumed in the Perkins model that the F layer is a slab of plasma 
in which the density and the derivative of the density go to zero at 
the endpoints. This is an oversimplification of the actual shape of 
the F layer, as will be shown later. 

The Perkins model is rederived now in the laboratory frame 
of reference including the boundary conditions of the F layer "slab" 
and also recombination. The other assumptions of the model are 
detailed in Chapter TI. 
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The equation of motion for the ions in the F region. is: 
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(3,1) 


The ion velocities parallel and perpendicular to B are obtained from 
Equation 3 . 1 . 
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In these equations, _V 0 is the electric field measured in the stationary 
frame of reference, and not the resultant electric field used by 
Perkins (1973). 

The perpendicular electron velocity is obtained in'a similar 
manner. Neglecting collisions and gravity, one ohtains: 
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The ion continuity equation is: 


Qli + V. A - nv^. + (nv z .) = -Pn 


a t 


(3.5) 


The quantity (3 is the loss rate due to recombination. Equation 3.5 
is integrated along a field .line, which goes down-in the northern 
hemisphere. 
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The assumption of plasma neutrality allows one to subtract the 
electron continuity equation from the ion continuity equation and 
integrate this result to obtain the following condition: 


*.l 




(3i7) 


It has also been assumed in obtaining Equation 3.7 that the electron- 
and ion velocities parallel to B are equal. This is a good assumption 
because high electron mobility keeps the electrons with the ions,' 
thereby preventing a space charge buildup along .field lines in -the 
F region. 

Equations 3.7, 3.6 and 3.4 are combined to give the equation 
describing the time evolution of the field integrated density. 
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The quantity N F in Equation 3.8 represents the total plasma content 
of a flux tube. It will change in time as a net flux of particles enters 
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or leaves a flux tube, as particles are lost by recombination, or as- 
the E x B drift causes a gradient to move into or out of the observed 
location. 

Equation 3.8 differs from that of the Perkins model (Equation 

2.10) in a number of ways. First, there is the inclusion of boundary 

terms which represents particle fluxes into or out of the flux- tubes. 

The Perkins model, by assuming n = 0 and ^ — = 0 at the end points,, 

assumes that there is no particle flux at the boundaries. Second, 

recombination is included. Third, the quantity V $ in the Perkins 

model is the electric field as seen>in a moving frame of reference, 

and implicitly contains the neutral wind. 

The fourth difference is the gravity term in Equation 2, 10. 

This term was carried through in Perkins' derivation because he did 

not write the continuity and current equations to the same order in 

v. /co .. If he had done that, he could then substitute the V • J = 0 
in ci — — 

equation in the more complete continuity equation, and the gravity 
term would vanish. This result must be obtained to be consistent 
with Equation 3.7. 

The Pedersen conductivity, which is the conductivity for currents 
parallel to _E and perpendicular to JB, is defined as: 
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In the F region, v.^<< w c j/ so tb 3 -*- cari be approximated by: 
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To obtain an equation describing the time evolution of the F region 


Pedersen conductivity. Equations 3.2 and 3=3 are substituted into 

v. 

Equation 3=5 and the resultiis multiplied by -g- — and integrated 

2 2 C ' - 

along a field line. Terms of order v. /co . are neglected. One 

m ci ° 

obtains: 
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where: 
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Because Equation 3. 10 was derived in the laboratory frame, it 
is different from Equation 2. 11 of the Perkins model. The trans- 
formation of Equation 2. 10 to the moving frame is accomplished by 
making the substitutions 
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where v 1 and V <p' are the velocities and electric fields in the moving 
frame,. Also, the. relation',. V ^ sin D = cos ^ has been used.' 

Equation 3.10 describes the time evolution of the field integrated 
Pedersen conductivity in the nighttime E region. The gradient term 
on the LHS represents a convective change in 2^, due to the x B 
and _E^x_B_ drifts of gradients in the Pedersen conductivity. The 
electric field and neutral wind terms express the. change in the 
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integrated Pedersen conductivity as the F layer rises and falls due 

to changes in these forces. The term containing NA represents 

-b 

the increase in Pedersen conductivity as the plasma falls down a 
field line due to gravity. The terms .(v 3?) represent the change -in 
conductivity as a net flux goes in or out of the flux tube. The 
physical significance of the boundary term containing {n^, - n^) 
is unclear. It has been labelled as a diffusion type term because it 
comes from the diffusion term in Equation 3.3. It probably represents 
the loss or gain of Pedersen conductivity as particles enter or leave 
the flux tube because of the vertical density gradient. The recombina- 
tion term represents the decrease in Pedersen conductivity as plasma 
particles are lost by recombining. 

3.2 Introduction to the Data 

The data used- in this thesis were taken at the Arecibo 

Observatory, Arecibo, Puerto Rico, using the incoherent-scatter 

technique (Evans, 1969; Behnke and Harper, 1973). Incoherent- 

scatter is a technique that measures the radar scattering from the 

ionospheric plasma. It is possible to obtain plasma drift velocities, 

ion temperatures, and density profiles using this technique. 

The Arecibo data are ideal for studies of mid- latitude F for 

s 

the following reasons. First, the Observatory provided high 
resolution data for this study. Second, the probability of observing 
F g at Arecibo is high, so it was possible to study nights with F . 

Third, the dip angle at Arecibo is 50°, which makes it ideal for 
studying mid- latitude models. On the other hand, it has the 
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disadvantage that models of the neutral wind can not be used at 
Arecibo because anamolous pressure behavior prevents a calculation 
of the east-west component of the neutral wind. The s ignif icance 
of this will become apparent in Chapter IV. 

This study concentrates on, four nights: September 17, 1974; 
November ^, 1974; May 18, 1975; and October 15, 1975. Ionograms 
for these nights, also provided. by Arecibo, were inspected. It was 
found that F g occurred on. the Setpember and May nights, but not on 
the November and October nights. Zinchenko (-1976) states that F 

s 

first appeared at 2150 on September-17, 1974 and at 1945 on May 18, 
1975. A personal inspection. of these ionograms- indicates that there 
is considerable uncertainty associated with trying to assign definite 
onset times,. This uncertainty is as much as several hours past the 
onset times given by Zinchenko (1976), and is primarily due. to. the 
traces being too weak to read. 

The data tape (provided by Dr. Robert Harper of the Arecibo 
Observatory) contains density profiles, with a density measurement 
every 1.2 km from 82 to 156 km, and every 6 km from 156 to 475 km. 
There are ion drift velocity profiles, with measurements every 29 km 
from 201 to 491 km, yielding 11 such- values, and a signal to noise 
(S/N) ratio associated with each-measurement. A S/N ratio of 
greater than 1.0 is required for the velocity measurements to be 
meaningful because the statistical error in these measurements when 
the S/N ratio is less than 1.0 increases rapidly with decreasing S/N 
ratios below 1.0 (Evans, 1969). 
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The plasma velocity measurement's available at each of the 
11 heights are v^g (positive north), (positive east) and v^ 

(positive with height). Conversion to the coordinate system given in 
Figure 2 is accomplished by the following transformations: 


v = v TT cos D - v, Tr , sin D 
x H NS 


V = V 

y E¥ 


v z = -v H sinD - v NS cos D 


The tape also contains ion temperature profiles from 137 to ■ 
492 km, with values every 29 km. 

The above data are given approximately every 20 minutes 
through the night, 

A computer program was also provided by the Arecibo 
Observatory which produces ion-neutral collision profiles using a 
neutral atmosphere model. Since these collision frequency profiles 
were used in this analysis, a brief description of the model is given 
in the Appendix. 

In this chapter, the comparison of the equilibrium values of 
Sj., and predicted by the model equations (Equations 3.8 and 
3.10) will be compared with the data described in the previous 
section. The model equations describe the behavior of N_, and S 

r r 

as changes in other ionospheric parameters occur. For example, 
Equation 3.8 states that the integrated density will change as more 
particles enter or leave a flux tube via the field lines, and as particles 
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are lost through recombination. Equation. 3. 10 states that the 
integrated Pedersen conductivity will change as the F layer rises 
and falls due to gravitational, collisional and electrostatic forces. 
Changes in the integrated conductivity will also occur as particles 
enter or leave a flux tube along the field lines, or are lost due to 
recombination. 

To completely analyze Equations 3.8 and 3.10, knowledge of 


the following quantities is needed: <B <p,>, N H, E , 2 

f .D r oy i: 

V NZ’ ^in^T’ ^in^B’ ^ V zi^T’ ^zi^B’ ^ — -*■" ^F’ and — -l- N F' Most 

of these can be inferred from the data as will be shown in the next 


section. 


Before starting, the analysis, it is useful to look at the density 
profiles shown in Figure 7, The F layer is well defined, but it 
is immediately seen that the assumption of Perkins (1973) that the 
density goes to zero s at the top and bottom of the layer, is an 
oversimplification. It is also seen that the shape of the layer does 
not change greatly during the night. The major shape changes occur 
at the bottomside of the layer, when the layer is low. This is 
attributed to recombination. Because the recombination rate depends 
on the neutral atmosphere density, the loss of particles due to 
recombination will be greatest at the lower altitudes, and particles 
will be removed from the bottom of the F layer at a faster rate than 
at the top of the F layer. This causes the increase in the density 
gradient from 19.82 to 22. 17 at the bottom of the layer. 

It is also noticed that the F layer moves up and down throughout 
the night. The motion of the layer is usually defined by the change 
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in time of the height at which the F region density, is. a maximum (h.m). 

According, to density profiles shown in Figure 7, hm is 320. km at 

19. 82, 380 km at 22. 17, and moves back down to 336-km at 1 . 30. 

According to Behnke and. Kohl (1974), the motion of ' the layer is governed 

by. the east-west electric field and the north-south neutral wind. This 

will be examined in the next section. 

It should. also be noted that the density profiles shown in Figure 

7 are for a night with F . The presence of the irregularities that 
^ s 

cause F do not seem to have an effect on the density profiles, 
s 

3.3 Analysis of the Data.- Fixed Limit Integration 

■The initial step in the data analysis must be a selection of the 
boundaries of the F region. • A method. whereby the boundaries of the 
F region are fixed regardless of the motion of the layer is presented 
in. this section. An alternate method is presented in the next section. 

The choice. of fixed boundaries is limited a great deal by the 
data. While there are density profiles and collision frequency profiles 
extending down to 82 km, the velocity data only extends down to 200 km. 
To make matters worse, the S/N ratio is dependent on the density of 
particles at the height of the measurement. When the layer rises, 
the density of particles at 200 km can be quite low as shown in 
Figure 7. The layer does rise to quite high altitudes, and can. remain 
high for long periods of time. For this reason, the S/N ratio for 
measurements taken at the two- lowes t altitudes (201.6 and 230.6 km) 
fall below 1.0 a significant amount of time. Thus, the third height, 
which is 259. 6 -km, is chosen as the bottom of the F layer. The 
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tenth height, which is 462.4 km, is chosen as the top because the 
density data only goes to 475 km. 

A recombination coefficient is calculated using the Stubbe 
model as referenced in Behnke and .Kohl (1974). 

P = 1.5x10 ^ n (O^) + 5.0 x 10 ^ n (N^) 

The oxygen and nitrogen densities are calculated using the models 
described in the Appendix. 

The integrals needed (N , 2 , <P>, P ) are calculated using 

_c i? 

Simpson's Rule. The range of integration was from 255. 0 to 
463.6 km; there are 36 points in this interval. 

Once the integrals are calculated, complete density and 
collision frequency profiles are no longer needed for the analysis of 
Equations 3.8 and 3. 10. Because it takes a large amount of computer 
time to simply read the data tape without doing any calculations 
(around 60 seconds), it was decided to take only the data needed 
from the tape and punch it on a disk system. This is a very 
inexpensive storage system that is easily accessible, in contrast 
to the tape which must be hand loaded for each run. 

The following data were stored, on disks; 1) hm, < P >, P 
2_, and T.„„ for each time, and 2) S/N ratio, VHA, YNS, n, 

0 ii ► 

, and q ■ • ■■ for each time at- each of the 10 heights. The -quantity 

hm is the height at which the density is a maximum, and T,™. is the 

bit 

exospheric temperature defined in the Appendix. The quantities 

S/N, VHA, VNS, n , v. , and 4^-are the S/N ratio, vertical ion 

e in a z 

drift velocity, nor th- south ion drift velocity, electron-density. 
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ion-neutral collision frequency, and electron density ‘gradient-'albng 
a field line, respectively, for each of the ten heights . 

A..slight problem -arises because- the density -measurements 
are not given at the same heights as the velocity measurements-. An 
example is shown below: 


Density Heights 
191.2 

197.0 
202.8 
208.6 

214.4 

220 . 2 

226.0 
231, 8 
237. 6 

243.4 


Velocity Heights 


201.6 


230.6 


To interpolate the density measurements to the same heights 
as the velocity measurements, a linear least squares fit was performed 
on the natural logarithm of the density. This is justified because 
particles in diffusive equilibrium will distribute themselves exponentially 
with altitude. Five points centered around the altitude at which an 
interpolated density was required were used in the fit. Thus, to 
interpolate the density at‘201.6 km, the fit was performed on the 
densities at 191.2, 19-7,0, 202.8, 208.6, and 214.4 km. The slope' 
of line is the derivative of the logaritm of the density, so the density 
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gradient is easily calculated. The same method of interpolation 
was used for the collision frequencies. 

September 17, 1974, which was the first day on the tape, will 
now be analyzed in detail using the quantities calculated. above. 

The motion of the F layer peak during the night of September 17 
is shown in Figure 8. It is seen that the peak moves up and down 
periodically, reaching a maximum of almost 400 km at 2230 and a 
minimum of 230 km at 0400, This behavior is typical of nights 
previously reported by Behnke and Kohl (1974). It becomes obvious 
at this early stage of the analysis that the choice of boundaries at 
259. 6 and’462.4 km is going to lead to trouble during those times 
when hm is near its minimum and its maximum. At 2230, when the 
layer is high,- the density of particles at 259.6 km will be so low 
that the S/N ratio will fall below 1.0. At 0400 when the peak has 
fallen to 230 km, the lower boundary is actually above the peak. 

This is not really in the spirit of the flux tube model as derived by 
Perkins (1973)'; one should at least include the maximum of the 
function that is being integrated if the assumption that n -*• 0 at both 
limits of integration is to have any meaning at all. The alternate 
scheme presented in the next section represents a solution to this 
problem. However, the fact remains that Equations 3,8 and 3.10 
should predict the behavior of N^, and 2^, for any limits of integration 
if the boundary conditions are properly, treated, 

Behnke and Kohl (1974) state that the motion of the peak is 
determined by the electric field and neutral winds. The Perkins 



hm (Km) 
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model also assumes that the nighttime F layer is supported by 
the eastward electric field and the neutral wind. 

The electric fields ar.e inferred from the ion drift velocity ' 
measurements. This can be explained with the help of Equation. 3 . 2, 
and the following approximate values for F region parameters: 


T 

in" 13 
10 ergs 

M = 

- 23 

2x10 “ g . 

CO . = 

CI 

2 x 10^ sec * 

g 

in 3 2 

10 cm - sec 

V, 

_ in 

CO . 
Cl 

-2 -3 

10 - 10 J 


It is clear that terms of order v. /<o . will not contribute much to the 

in ci ... 

perpendicular ion drift velocity. The _g x z term will only contribute 
a term on the order of .05 m/s to the velocity. The gradient term 

7 

is of the order — ^ cm/sec, where H is the characteristic 

h p p 

length of the density gradient perpendicular- to the field lines. This 
length would have to be on the order of 1 km for the density gradient 
term to even contribute 5' m/s to the perpendicular drift velocity. 

In fact, the characteristic length of the density gradient in the F 
region is in the range 50 to 100 km and the drift velocities are 
generally much larger than 5 m/ s. One can conclude that the 
dominant cause of perpendicular ion drift is the electric field: 


— Di 

c „ A 

•g- Exz 

V Dx = 

X E oy . 

V Dy = 

- 1- E o* 
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It should be noted that electric fields in the F region are usually- 
inferred by this method. It is both conventional and convenient to ■ 
represent the magnitude of ele.ctric fields in terms of the drift 
velocity they would cause rather than the actual unit of statvolts/cm. 
Thus, electric fields are multiplied by c/B and are given in units - 
of meters/sec . 

It was stated that there are drift velocity measurements for 
10 heights with a S/N ratio associated with each measurement. An 
average F region drift velocity is obtained by multiplying each of 
the 10 drift velocities by the square root of the S/N ratio associated 
with each measurement, summing the products, and dividing by 
the sum of the square roots of the S/N ratios: 


10 



i= 1 


A measurement is thrown out if the S/N ratio is less than 1.0. In 

this manner,- average F region values are found for v and v n , and 

ux u y 

the equilibrium electric fields inferred from the averages. 

The eastward electric field during the night of September 17, 
1974 is shown in Figure 9. There is some correlation between the 


direction and magnitude of the eastward electric field and the motion 


of the peak shown in Figure 8, but it is not good enough to say that 
the electric field is the dominant lifting force. This fact is especially 



17 , 1974 
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clear if one looks at the times before 2200 when the layer Is rising 
rapidly. The electric field is in the westward direction, meaning 
it is causing ion motion downward, 

The north- south neutral wind, which can literally blow ioniza- 
tion up the field lines, is the missing factor. It will now be shown 
that the neutral wind and the electric field are the dominant factors 
governing the 'motion, of the F layer peak, further substantiating the 
work of Behnke and Kohl (1974). 

The neutral wind. must also be inferred from the drift velocity 
measurements. Equation 3.3 provides a direct' solution of the neutral 
wind along the magnetic field (V^.^.) with v^. being the parallel ion 
drift velocity given by the data, 

v _ g sin D 2T 9 n 

NZ V zi v. Mnv. 5 z 

in in 

i 

The neutral wind is assumed to have no vertical component. 

. Using the available data, one can solve for a value of Y AT7 at 
each of the 10 heights. Since an average F region value, is desired, 
a logical first step would be to solve for a value of V at each 

IN /Lt 

height, and then find an average by weighting with the square root 
of the S/N ratio as described previously. This method does not 

. r\ 

yield good results, however, -because small fluctuations in • 

9 z 

at high altitudes are greatly magnified by 1 / v n * 

An alternative method used by Behnke and Kohl (1974) is to 
evaluate V at the peak, where 9 n/ 9 z- is zero and the S/N ratio is 

liZ 

large, and assume that this is the average velocity of the neutral wind. 
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A comparison of the results of the two methods .is shown in Figure 10. 
The general behavior is the same, but the value of V.^^; obtained at 
the peak shows far less fluctuation. For this reason, the use of 
at hm will be adopted. 

The neutral wind is added to the electric field to produce "the 

resultant electric field. The value of the eastward component of the . 

resultant electric field (E——^ ) is shown in Figure 11; hm is plotted 

-KlLo Y 

on the same graph to show the excellent correlation between hm and 
^RESY' There are only two periods in which the motion of peak 
can not be explained by the resultant eastward electric field. These 
are the time periods 1900 and 1930 and 0400 to 0430. In these two 
periods the layer is rising, although the field is westward. Even 
though-the magnitude of the westward field is decreasing steadily 
until the field becomes eastward, a westward field of any magnitude 
should produce a downward motion of the peak, reinforcing the effect 
of gravity. The behavior of the peak in these two periods can be 
explained with the help of Figure 7. It is seen that recombination has 
removed more particles below the peak than above the peak between 
the times 19.82 and 22. 17 shown in Figure 7. This will cause the 
peak to rise, even in the absence of bulk plasma motion upward. The 
same conclusion holds for the time period after 0400, The layer is 
quite low, so recombination would be expected to be important. One 


can thus conclude that the addition of recombination to the model is 


not a trivial extension, because it is required to explain the peak 
behavior during certain times of the night. 



(m *«£' ) 



Figure 10: Time variation of the neutral -wind at peak and the neutraTwind 
weighted by the square root of the S/N ratio. 
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'Figure 11: Time variation of the resultant eastward electric 
field and hm for September 17, 1974. 
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Before continuing- with the analysis of the night of September 17, 

it is useful to show that the s.ame correlation-betwe,en an< ^ h™- 

is -found on the other nights (Figures -1-2, 1-3 and- 14). -It -is seen 

that the resultant electric field on, the November night, which, is a 

night without F , is much smaller in magnitude than the field for the 

two nights with F (May 18 and September 17). It also shows far less 
s 

time variation. The field during the October night, which is also a . 

night without F , shows smoother time variations than the two nights 
s 

with F , but its behavior is not as smooth as the November night, 
s 

One can not draw any conclusions based on this. However, it should 

be noted. that the resultant eastward electric field provides a stabilizing 

force for the Perkins model of F ; it is the northward field that causes 

s 

the instability to develop. The northward electric fields for the four 
nights will be examined at a later time. 

The analysis so far has shown that the correlation-, between the 
eastward resultant electric field and the motion of the F" layer peak 
is good, with the exception of times when the- layer" is low and 
recombination must be added. One could develop a model describing 
the motion of the peak by including only the eastward electric field, 
the neutral wind, and recombination, as essentially was done by 
Behnke and Kohl (1974). However, the Perkins model is a flux tube 
model, whereas a model describing the motion of the peak would be 
primarily local. 

The Perkins model states that the field integrated Pedersen 
conductivity (2^) will change in time as the layer is moved up and 
down, by gravitational, electrostatic, and collisional (neutral wind) 
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Figure 12: Time variation of the resultant eastward electric 

field and hm for November 9, 1974. 
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Figure 14: Time variation of the resultant eastward electric 

field and hm for October 14, 1975. 
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forces. If the layer shape remains reasonably constant, as it does 

(Figure 7), then the time behavior of 2^, should be anti- correlated 

with the behavior of the peak, because v. / w . decreases with an 

. in c.i 

increase in altitude. As the layer goes down, the integrated conducti- 
vity should increase because it depends on /w^. 

It is seen by looking at Figure 15 that 2^, and hmare indeed 
anti- correlated. When the layer rises, 2^, decreases, and when it 
falls 2^, increases. The time period that represents an exception 
to the above statement is from 0230 to 0400, when the layer has fallen 
to 230 km but 2^, has not increased. This, is primarily due to choosing 
the lower limit of integration at 259. 0 km. When the peak falls 
below this value, as it does at 0400, a sizable number of particles are 
not counted in the integration process. The value of 2^, also depends 

on N_, the total number of particles in the flux tube. Thus, if N 
_b F 

decreases, 2^., will also decrease even if the layer is not moving up. 

Failures in. the anti- correlation between 2„ and h can usually be 

F m 3 

attributed to changes in N . The. alternate scheme to be presented of 
^f an d hm for the other nights is shown in Figures 16, 17. and 18. 

The anti- correlation for the May night is the weakest. 

It has now been shown that the resultant electric field is the 
dominant mechanism governing the height of the peak, and that the 
height of the peak is anti- correlated with the value of the integrated 
Pedersen conductivity. It can be concluded that the equilibrium 
solution of the Perkins model, in which the value of the field integrated 
conductivity depends on the resultant eastward electric field and the 
integrated density (Equation 2.11) is physically meaningful. 
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Figure 15: Time variation of the F region integrated Pedersen 
conductivity and hm for September 17, 1974. 


Jun (km) 


(xl0 8 cm sec 1 ) 


NOVEMBER 9,1974 


2000 



“1 1 1 1 1 r i 1 1 1 1 1 1 r 

18 20 22 24 2 4 6 

LOCAL TIME (HOURS) 

Figure 16: Time variation of the F region integrated Pedersen 
conductivity and hm for November 9, 1974. 
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Figure 17; Time variation of the F region integrated Pedersen 
conductivity and hm for May 18, 1975. 
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Figure 18: Time variation of the F region integrated Pedersen 
conductivity and hm for October 14, 1975, 
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Equation 3. 10 has a lot of extra terms in it that would, not be expected 
to be important based on the above conclusion, at least for most of 
the night. One exception is recombination, which has been found to 
be important when the layer is low. The effect of the extra terms 
will be discussed. laterin this chapter. 

The behavior of the integrated plasma density' (Np), is shown 

in Figure 19, It is seen that the integrated density decreases steadily 

until around midnight, when, it increases. It then holds constant for a 

few hours until it starts to decrease again after 0200, Since 

recombination will cause a steady depletion throughout the night, the 

increase in must come either from a density gradient drifting 

into the region of observation, or from a flux of particles coming in 

at the boundaries of the flux tube. The decrease of through most 

of the night can be explained by recombination, but the midnight 

increase can only be explained by more particles coming' into a flux 

tube than are lost by recombination. It is the assumption at this 

stage that a net input of particles at the boundaries of order 3 x 10^ 
-2 

cm over a period of 40 minutes is the dominant mechanism that 
accounts for the input of particles. This requires an input flux of 

„ g ^ 2 _ J 

roughly 2x10 cm - sec . The question of the density gradient 
will be considered later. 

The fluxes are determined at each of the 10 heights simply by 
multiplying the bulk ion drift velocity along the magnetic field by 
the density at the height in question. If the S/N ratio is less than 
1. 0 at any height, the flux is set equal to zero for the purposes of 
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plotting. The fluxes at each of the 10 heights for the night of 
September 17 are shown in Figures 20 a-j, A positive flux is 
downward along the magnetic field. 

It is immediately seen that the fluxes at the lower heights are 
generally smaller than at the higher heights. The velocity that 
contributes to the flux is a combination of the diffusion velocity and 
the neutral wind velocity (cf. Equation 3,3), At the higher altitudes, 
1/ v.^ is quite large so. the diffusion velocity is expected to be large. 
On the other hand, n is decreasing with the result that the flux at tfe 
higher altitudes should not change significantly from one altitude to 

e 

the next since recombination is negligible. In fact, the fluxes 
calculated by these data do change over adjacent heights. This is 
illustrated in Figure 20 k. Also, the fluxes vary a great deal in time. 
These two facts raise serious questions about the accuracy of the 
fluxes. Further indications of the inaccuracy of the fluxes will 
appear when. the fluxes are used in Equations 3,8 and 3, 10, 

Because of the large amount of time and space variation of the 
fluxes, an attempt was made to calculate the fluxes directly from 
Equation 3.3. This method does not work, however, because of 
the 1/v dependence of the diffusion velocity. Even with a least 
squares smoothing of density data near the top, there are still small 
fluctuations in 3n/ 3z that are magnified by the l/v\ factor. The 
resulting fluxes calculated directly by Equation 3.3 showed even 
larger variations in both time and space than those calculated with 
velocity data. 
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Figure 20c: Time variation of the flux at 259.6 km. 
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Figure 20e: Time variation of the flux at 3jj;7. 6 km. 
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Figure 2 0 j : Time variation of -the flux at 462.4 km. 
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The quantity desired in Equations 3.8 and 3. 10 is the net flux 
into a flux tube, i.e. $ T - i> B< When the S/N ratio is below 1. 0 for 
measur em ents at 259.6 km , the velocity at the next higher height 
that has a S/N ratio of greater than 1. 0 is used with the density at 
259.6 km to calculate the flux. The net flux into a flux tube is shown 
in Figure 21 . 

It is obvious that the net flux shown in Figure 21 will not 
explain the behavior of the integrated density. The period just 
before midnight, when the density is increasing, is one when there 
is a large amount of particles leaving the flux tube via the field lines. 

8 - 2 _ i 

A flux into the flux tube of 2 x 10 cm - sec is required to explain 

8 

the density increase; the data shows , a flux of greater than 5x10 

leaving the flux tube. The net out- flux of particles in the period 

12 - 2 

from 2300 to 0100 represents a loss of around 2 x 10 particles cm 

This would certainly be seen in Figure 19 if such an outflow did 
* 

indeed exist. 

Because the fluxes calculated at the boundaries (259. 6 and 
462.4 km) were so different, causing a much larger net flux than the 
integrated content data indicated, it was initially thought that narrowing 
the integration interval would alleviate the problem somewhat. The 
behavior of the fluxes at adjacent heights (Figure 20k) indicates that 
this procedure is not going to help much. There are other reasons 
for choosing different boundaries, however, and they will be discussed 


in the next section. 
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Figure 20k: Time variation of the flux at 404.5 km and 433. 5 km. 
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3.4 Analysis of the Data --Peak-Centered Integration 

The problems associated with the fixed limit integration scheme 
are two-fold. The first is. that when the peak is at a maximum or 
minimum height, the S/N ratio at one of the boundaries may be less 
than 1.0. The second is that the peak may actually fall below one 
of the boundaries and the integration will not include a sizable portion 
of the layer. The method of analysis presented in this section eliminates 
both of these problems by narrowing the interval of integration and 
requiring that the interval follow the peak up and down. By narrowing 
the interval (to 90 km) and centering it around the peak, it is insured 
that the S/N ratio will be high at the boundaries. In addition, it was 
hoped that perhaps the net flux would show better agreement with the 
data on N„. 

The net flux calculated by the peak— centered method is shown 
in Figure 22. As is seen, it still shows a great deal of time 
variation, and the net outward flux before midnight is still there. 
Consequently, this problem is still unresolved. 

Now that a method of analysis has been chosen which is -in the 
spirit of the Perkins flux tube model. Equations 3.8 and 3. 10 can 
be compared to the data. Assuming no horizontal gradients in N or 2, 
the two model equations are of the form: 


dx 

dt 


-A(fc) x + S(t) 


where A(t) and S(t) are determined from the data. The solution to 
this equation is 
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The integrals in this equation are calculated, using the trapazoidal rule. 

The solution of Equation 3.8 is' compared with the actual values 
of N_(t) calcxilated by the data in Figures 23 and 24. In Figure 2-3’ 
the fluxes are not included, while in Figure 24 they are included. 

It is immediately seen that the model with only recombination added 
represents, the data quite well. When the net flux is added, the model 
breaks down completely. The fluxes are just too large and turbulent 
compared to what the data for N_ predicts they should be. The fact 
that the- model with only recombination added predicts the density 
behavior fairly accurately suggests that the measured fluxes are 
much too high. The' same conclus ionsis reached when the other 
nights are analyzed (Figures 27, 30 and 33), although the agreement 
between the model with recombination added is not as good as it was 
during the September night (Figures 26, 29 and 32). This indicates 
that there is some flux entering or leaving, but it is not the flux 
given by the data. 

The flux terms in the equation describing the behavior of 2 
* > * - _b 

(Equation 3.10) also give trouble, as will now be demonstrated.' 

Terms -were successively added to the model by firs t' including only 
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Figure 23: Time variation of N^, obtained frc 
Equation 3.8 including recombina 
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Figure 24: Time variation of obtained from the data, and the solution of 
Equation 3.8 including recombinatLqn.-and'the flux. ■ 
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Figure 26: Time variation of N^, obtained from tHe data', and the solution of 
Equation 3.8 including recombination only.. 
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Figure 27: Time variation of N_ obtained from the data, and the solution of 
Equation 3.8 including recombination and the flux. 
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Figure 29: Time variation of Np obtained from the data, and the solution of 
Equation 3.8 including recombination only. 
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Figure 30: Time variation of TsL obtained from the data, and the solution of 
Equation 3.8 including recombination and, the flux. .. 
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Figure 32: Time variation of N^, obtained from the data, and the solution of 

Equation 3.8 including recombination only. ' 
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Figure 33: Time variation of Np obtained from the data, and - the solution of 
■ Equation 3.8 including recombination and the- flux. 
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the electric field, neutral wind, and gravity terms, then adding 
recombination, then the diffusion term, and. finally the flux terms. 

The results for September 17, 1974 are shown in Figures 34-37 and 
compared with the actual values of calculated by the data. It 
is noted that the flat regions merely mean the solution is off the 
chosen scale for the plot. It is evident that recombination is important 
when the layer is low, and that the addition of the diffusion term 
helps the model slightly. The addition of the fluxes to the model 
again causes it to break down, except right at the beginning. 

The other nights show similar behavior, as shown in Figures 
38-43. The one exception is the May night, when the addition of the 
fluxes did not hurt the model greatly. 

The issue of why the flux data, which have been assumed by 
most researchers to be accurate at Arecibo, do not agree with the 
rest of the data is a very important one. It will be examined in the 
next section. 

The Perkins model of F g assumes a quas i- equilibrium in 
which the time scale of the change in the equilibrium integrated 
Pedersen conductivity is much longer than the time scale of the 
development of F . As has been shown, the equilibrium- value of 
the field integrated Pedersen conductivity changes rapidly through 
the night. To test the assumption that there is a quasi— equilibrium, 
the Perkins equilibrium condition (Equation 2. 12) is used to obtain 
the equilibrium value of Z^. The results are shown in Figures 44-47 
where the "model" solution includes only 


This is consistent 



(x I0 8 cm sec 1 ) 




SEPTEMBER 17,1974 
E« s y ONLY 



LOCAL TIME (HOURS) 






X F (x 10® cm sec 1 ) 


SEPTEMBER 17,1974 
DIFFUSION ADDED 



100 




2 f (x 10® cm sec 1 ) 


NOVEMBER 9,1974 
DIFFUSION ADDED 


2000-1 


I500H 


1000 


500 -f 



0 *+* 

. 18 


T 


“i 1 1 1 1 — — j — ■ — i 1 r 

20 22 24 2 4 

LOCAL TIME (HOURS) 


t r 

6 


Figure 38: 


Time variation of 2 obtained from the data, and the solution of 

Equation 3.10 inciucung E-.^-,, gravity, recombination and' 

ivJtLiO 

diffusion. 
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Figure 39: Time variation of 2 obtained from the data, and the solution of 
Equation 3. 10 including E p gravity, recombination, diffusion 
and the flux. K 
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Figure 40: Time variation of obtained from the data, and the solution of 


Equation 3.10 inclucfing 
diffusion. 
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Figure 41: Time variation of 2 obtained from the data, and the solution of 
Equation 3.10 including E RTr c> gravity, . recombination, diffusion 
and the flux. i ^ 
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Figure 42- Time variation of S obtained from the data, and.,the solution of 
Equation 3.10 including E RJr _, gravity, recombination and 
diffus ion. 
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Figure 44: Time variation of 2 obtained from the data, the solution of Equation 3. 10 
including E RT?C , and gravity, and the quasi- equilibrium, solution of the 
Perkins model (Equation 2. 12), 
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Figure 45: Time variation of obtained from the data, the solution of 

Equation 3.10 including and gravity, ..and the quasi- equilibrium 

solution of the Perkins model (Equation 2. 12). 
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Figure 47: Time variation of S obtained from the data, the solution of 

Equation 3. 10 including E „„ and gravity, and the quas i- equilibrium 
solution of the Perkins model (Equation 2. 12). 


Ill 


112 


with Perkins. It is seen that the agreement is reasonably good for 
the September, May and October nights, especially when the layer 
is high. The qua si- equilibrium predicted by the Perkins model is 
not good at all for the November night. Curiously, this is a night 
without F g and with' rather smooth behavior of the equilibrium 
parameters . 

The assumption of a quasi- equilibrium is really a comparison 
of time scales: 


.o 


9 2 

c 

9 t 


<< 


■where y is the growth rate of F . Around the onset time for the 

3 

September night (2150), 


9 2 

c 

9 t 


3600 


sec 


-1 


The time scale for the development of F is assumed to be from 
10 minutes to 30 minutes, so the assumption of a quasi- equilibrium 
is. indeed a good one. For the May night, 2 is decreasing quite 
rapidly around onset. The quasi- equilibrium predicted by the Perkins 
model is not close to the actual value. - 

It is the conclusion of this analysis that great care must be 
taken when a quasi- equilibrium of the form in the Perkins model is 
chosen. It is a fairly good assumption on some nights, but- on others 


it is not. 
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3.5 Fluxes and Gradients 

A great deal has been- written so far about, the disagreement 
between the flux data and the rest of the data. The degree of 
disagreement will be demonstrated in this section, and some alternate 
explanations proposed. 

One can solve for the flux required to give agreement with 
the rest of the data by solving Equations 3,8 and 3. 10 simultaneously, 
in the following manner: 


§? - a = A 

T B • 


(3.12) 


(v in ) T $ T ^in^B . C 


(3. 13) 


The quantities A and C are found from Equations 3., 8 and 3-. 10: 


3 N, 


A = ~^r~ + < P > N 


C = 


co . B 
ci 

ec 


8Z f 

£_ q Sl- 

3 t B 


-22J-2- E s 
H oy F.- 


ec g sin D 


co . B 
ci 


N_ 


H - ■ ■ F 


- sinD .. Y 2 

H NZ F 


ec 2T sin D 


co . B M H, 

ci 


(n T - n B ) + ; P , Sj 


The simultaneous solution of Equations (3.12) and (3,13) yields the 
required values of and c& : 

1 JD 
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Aw - r. 



The results of this calculation, compared to the flux given by 
the data are shown; in Figures 48-51 for two of the nights. It is 
immediately seen that the; flux required is not even close to the flux 
given by the data. This result is - hardly surprising based on previous 
attempts to use- the. flux data. 

Much more, informative are the results of the calculations of 

the net flux required to satisfy Equations- 3. 8 and 3.10 compared to 

the flux given by the data. This:is shown- in Figures 52-55. It is 

seen that the flux. required to satisfy the N_ and 2 data are generally 

i* F 

much lower than. those given, by the data. In<addition, the required net 
flux. does not exhibit the large amount of time variation that is seen in 
the flux calculated by the data. 

An error analysis was performed on the flux .data. It is seen 
from Figure 52 that the flux required and. the flux from the data 
differ by a factor of 10 or greater for a large portion of the night. 
Assuming that the net flux required is on the order of 5 x 10 , while 

Q 

the net flux from the data is on, the order of 5 x 10 , one would have 

O 

to subtract an error of 4. 5 x 10 to obtain agreement. The error 
analysis proceeds in. the 1 following, way: 
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Figure 48: Time variation of top fiiuc, and flux required to satisfy 
Equations 3.8 and 3.10. 
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SEPTEMBER 17, 1974 



LOCAL TIME (HOURS) 


Figure 49: Time variation of bottom flux, and flux reqjiired to satisfy 
Equations 3.8 and 3. 10. 
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Figure 50: Time variation of top flux, and flux required to satisfy 
Equa tions 3 . 8 and 3.10. 
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Figure 51: Time variation of bottom flux, and flux required to satisfy 
Equations 3.8 and 3.10. 
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Figure 52: Time variation of net flux, and flux required to satisfy 
Equations 3.8 and 3.10. 
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NOVEMBER 9,1974 
FLUX DATA. 



Figure 53: Time variation of net flux, and flux required to satisfy 
Equations 3. 8 and 3.10. 
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MAY 18, 1975 



LOCAL TIME (HOURS) 

Figure 54: Time variation of net flux, and flax. required to satisfy 
Equations 3.8 and 3.10. 
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Figure 55: Time 1 
Equati 




12 3 ' 


A (NET FLUX 1 ) =\TT~ A (nv ) 

A (n v ) = n A v 

A{NETFLUX]= 'vFZ n Av 
n p* 1 x 10 5 

Ay= AMETJLgX L 4.5 x10 s = 32 m/sec 
<s /~2 n \l 2 1x10 

It has been assumed that the error is identical for both nor th- south' 
and vertical components of the velocity, and that errors in the density 
are significant compared to those in the velocities. If the uncertainty- 
in the vertical component is reduced to 5 m/sec, -the uncertainty in 
V NS wou ^ increase to 50 m/sec. 

It is difficult to come to the conclusion that the fluxes measured 
by such a high resolution instrument as the Arecibo radar are so 
erroneous. Before coming to this conclusion, -the possibility of having 
a density 'gradient produce an effect that would explain the 'behavior 
of Ey and N^, will be examined. The largest flux that would be 

g 

required to satisfy Equations 3.8 and 3. 10 is about 2.5-xTO . 

Assuming an electric field drift velocity of the order of 25 m/sec,' one' 
can solve for the density gradient required to yield a value’ lof 

O _ O _ 1 

2,5 x 10 cm - sec (cf. Equations 3,8). 
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0 , in 3 -1 AN 

2.5x10 cm - sec x 


Ax' 


9 r in 8 ' -2 -1 

= 2.5 x 10 cm - sec 


AN' , . ..5 -3 

1x10 cm 


Ax'. 


This has more meaning if written in terms of the characteristic 
length of the gradient: 


N 

AN/ Ax 


i in 12 " 2 
1x10 cm 

5 -3 

1 x 10 cm 


100 km 


Thus, if there were a density gradient of this magnitude (which is 

on the order of magnitude of the vertical density gradient) drifting 

horizontally at a speed of 25 m/sec, it could explain the behavior of 

Np. In the data examined, there was no evidence of a horizontal 

gradient of this magnitude. However, other investigators have found 

evidence of substantial gradients (Matthews and Harper, 1972). 

The above analysis merely states that a density gradient could 

explain the behavior of Np and Sp, assuming there was no flux. If 

one assumes. that the flux data is correct, then the effect of the 

density gradient Would have to counteract the effect of the flux and 

produce the required curves shown in Figures 52-55. This is -highly 

unlikely because some of the spikes seen in Figure 52 reach almost 
9 - 2 - 1 

1x10 cm - sec . To counteract these, and produce a resultant 

7 

flow of particles of around 1x10 would require a much larger density 
gradient, as is shown below: 
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y AN 
NET D Ax 


= 1 x 10 


AN 

Ax 


10 7 

2. 5 x 1CT 


-2 -1 
cm - sec 


cm - sec 


-1 


AN 

Ax 


.a 1 x 10 12 
4 x 10 5 


-2 
cm ' 

^3 


cm 


25 km 


This Is a very large density gradient. It is highly unlikely that 
horizontal density gradients of sufficient. magnitude would drift - 
horizontally in exactly the right direction to counteract the enormous 
fluctuations that are seen in the flux to produce the source terms • 
required as shown In- Figures 52-55. It is also highly unlikely that . . 
particles would be building up in the flux- tube, causing the horizontal . 
gradient, and then drift across the field lines. It is much easier for 
them to go straight on through the flux tubes along the. field, lines, thus 
causing no buildup. - . 

It is unfortunate that nothing more can be said about the 
density gradients or the flux data. There; is no direct data presently 
available for the horizontal density gradients, and it will probably 
be difficult to acquire such data. At the present time, the issue of 
the flux behavior remains unresolved. 

The tentative conclusion is reached that the flux data given 
by the Arecibo data have much larger uncertainties in them than has 
previously been assumed. 
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A possible explanation of the large errors associated with 
the flux is found in the recent work of Vickrey et al. , (1976). The 
reduction- of -the -raw A-r-ecibo -data. involves the assumption that -atomic 
oxygen is the single species of ions in the upper F region. Vickrey 
et al., (1976) have shown that there is a large amount of hydrogen 
ions in the F region. The large mass difference between hydrogen 
and oxygen could produce errors in the drift velocities calculated 
from the Arecibo data on the order of 80 m/s. 

3.6 Equilibrium Behavior of F Nights vs. Non F Nights 

s s 

The flux is somewhat more turbulent on F nights. (September 

s 

17, 1974 and May 18,' 1975) than on nights without F (November 9, 

s 

1974 and October 14, 1975) as is shown in Figures 22, 25, 28 and 31. 

It is also seen that the electric fields are more turbulent (Figures 9, 

56-62) on F’ nights, in agreement with Kelley (1972). 

There is no difference in the other equilibrium parameters 

such as the motion of the peak, 2^, or N^, during non-F g nights. 

The next chapter will examine the growth rate predicted by 

the Perkins model of E , and compare it to the data. 

s 
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Figure 57: Time variation of eastward electric field 
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Figure 58: Time variation of northward electric field. 
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Figure 59: Time variation of eastward electric field. 
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Figure 60: Time variation of northward electric field; 
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CHAPTER IV 


GROWTH RATE ANALYSIS 


In the last chapter, it was shown in Figures 44 and 46 that the 

equilibrium at the onset of F g predicted by the Perkins model was in 

agreement with the data for one of the days with F (September 17, 

s 

1974) but not the other day (May 18, 1975). In this chapter, linear 
growth rates predicted by the Perkins model and the Scannapieco 
model will be examined. Using the* Arecibo data, the conditions required 
to yield reasonable growth times will be obtained. In addition, the 
Perkins model will be extended to include the effects of short 
wavelength damping, recombination and E region coupling. 

4. 1 Comparison With the Data 

The linear growth rate predicted by the Perkins model is 


c E RES cos D 
B H 


sin 0/2 


where the direction of the propagation of the irregularities has been 
assumed to be that which gives maximum growth, i.e., a = 0/2. 
The linear growth rate of the Scannapieco model is 
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where k = k x is assumed in this model. 

— x 

Ideally, one would like to use the data to 'calculate growth 

rates predicted by the two models and compare the predicted rates 

during the onset of F with the predicted rates when no F occurred. 

s s 

If either or both of the models could predict the occurrence of F , 
long growth times or damping should be seen before onset, and short 
growth times at onset. 

This is not possible to do because of the lack of data on 

several critical parameters. One of these is the direction of 

propagation of the field aligned irregularities that cause F . This is 

, . . , . s . 

important because the Perkins model predicts damping if the direction 
of propagation is more northward than the direction of the equilibrium 
resultant electric field. This issue can be bypassed if one assumes 
that the propagation direction will be that which yields maximum 
growth, i. e. a = 9/2 for the Perkins model and a = ir/2 for the 
Scannapieco model. 

A much more serious problem is the lack of data on- the 
ea’st-west component of the neutral wind (V^^) at Arecibo, This is 
required to calculate the resultant electric field in the x direction, 
which is the very parameter that is required for positive growth in 
both models. It should also be noted that there are no models that 


accurately predict at Arecibo (Behnke and Kohl, 197.4). Finally, 

i 

the lack of information on 9 I/9y means that there are two. critical- 


parameters in the Scannapieco model which are unknown. Since 
there is only one unknown in the expression for the Perkins growth 
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expression for V^y, using the Arecibo data. This has been done for 

times around onset, and the results are shown in Tables 1-4. 

It is immediately seen by looking, at Tables. 1-4 that -the 

neutral wind required for growth times of 10 minutes is on the order 

of 200-500 m/s, and for growth times of 30 minutes, up to- 200. m/s. 

These are not unphysical results, as can be seen in Figures 10, 63, 

64 and 65 which show the neutral wind along the field line. 

The next point that is noticed in the October data- is that there 

is no difference in the neutral winds required for F g onset between 

days with F and those without F . However, the actual magnitude 
s s 

* •* r 

of the east- west neutral wind does seem to be different for nights 

with and without F g as seen in Figures 10, 63, 64 and 65. The winds 

reach much- higher velocities, at least along the field lines, on the 

days with F ,. However, on the May night, the wind reaches its peak 
s 

velocity at a time well beyond the onset time. 

The last point is that there is no great change in the neutral 

wind required between onset and the times before onset. This 

suggests that the other equilibrium parameters do not trigger F , 

. ' . s . 

and the quantity V^y must change before onset if the Perkins model 
is to predict the occurrence of F g . 

The conclusion can not be emphasized too strongly. Without 
an independent measurement of V^y, it is impossible to calculate 
growth rates predicted by the Perkins model, or the Reid or 
Scannapieco models, for that matter. It is not possible to infer that 

i > i 

because the values of V NZ are higher on F g nights that the values of 
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Table 1 

September 

17, 1974 




t = 10 Minutes 

t = 30 Minutes 

Time 

V 

NY 

|l ^o]- 1 

V 

NY 

(1 9 

i 2 a y 

S 9 


{ml s ) 

(km) 

(m/s) 

. ° (Ion) 

20. 90 

300 

-83 

150 

-50 

21. 15 

310 

-83 

160 

-53 

21.38 

330 

-83 

170 

-56 

21, 68 

320 

-77 

170 

-50 

ONSET AT 

21. 83 




21. 93 

380 

-63 

210 

-40 

22, 17 

400 

-62 

230 

-43 

22.47 

490 

-53 

300 

-31 

22. 72 

430 

-72 

260 

-43 

22. 95 

420 

-53 

260 

-32 

23. 25 

210 

-53 

48 

-31 

23. 50 

210 

-63 

42 

-38 

23. 73 

67 

-59 

-79 

-37 

. 03- ■ 

150 

-63 

9.2 

-37 


Table 2 

November 9, 

1974 



20. 00 

240 

- 91 

79 

-77 

20. 25 

240 

-111 

74 

-77 

20. 55 

210 

-103 

58 

-62 

20. 78 

210 

-114 

54 

-83 

21. 03 

210 

-100 

56 

-63 

21.33 

210 

-100 

58 

-62 

21.57 

220 

-100 

65 

-62 

21. 82 

210 

-100 

55 

-62 

22. 1-2 

210 

-100 

53 

-62 

22.35 

200 

-100 

53 

-63 

22. 60 

250 

- 83 

93 

-50 

22. 90 

260 

- 83 

110 

-53 
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Table 3 

May 18, 

1975 




fc = 

10 Minutes 

t 

: = 30 Minutes 



/ . 82 |-1 


i as 



f 1 ° 1 


. [ 1 o 

Time 

V 

NY 

s 3y 

V 

NY 

2 8y 


(m/s) 

1 ° (km) 

(m/s) 

' ° (km) 

19, 07 

260 

-143 

92 

-132 

19.32 

280 

-123 

120 

-100 

ONSET 

19. 50 




19. 63 

250 

-108 

110 

- 83 

19. 90 

250 

-103 

110 

- 77 

20. 15 

240 

-100 

99 

- 71 . 

20.47 

240 

-100 

98 

- 62' 

20. 73 

240 

- 91 

100 

- 62 

20. 98 

270 

- 83 

130 

- 59 

21.30 

250 

- 77 

120 

- 53 

21. 57 

250 

- 77 

120 

- 53 

21. 82 

220 

- 72 

84 

- 48 

22. 13 

230 

- 72 

96 

- 48 

22. 38 

220 

- 72 

85 

- 46 


Table 4 

October 14, 

1975 



18. 88 

230 

-100 

77 

- 62 

19. 13 

300 

-133 

130 

-106 

19.47 

260 

-118 

100 

- 91 

19. 77 

250 

-115 

99 

- 83 

20. 02 

240 

-108 

96 

- 77 

20. 35 

250 

-122 

100 

-111 

20. 62 

290 

- 77 

140 

- 50 

20. 87 

330 

- 83- 

' 170 

- 50 

21. 20 

340 

-132 

190 

-127 

21.45 

330 

- 72 

180 

- 4'6 

21. 72 

260 

- 56 

120 

- 33 

22. 03 

210 

- 72 

79 

- 48 
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Figure 63: Time variation of neutral wind. 
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- Figure, 65: Time variation of neutral wind 
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Y ny will also be higher and hence reasonable growth times. will be 
obtained. It is only possible to state that given the right values of 
Y ny -before onset and at onset, the Perkins model can predict the 
development of F g with reasonable growth times, and that the values 
required for Vj^ Y are not too large to be unphysical. 

The other two columns in Tables-1-4 show the characteristic 
size of the east-west equilibrium gradient (1/2^ / 3 2^/3 y) ^ that 
would be required to yield reasonable growth times-, in the Scannapieco 
model, assuming the values of Vj^ Y obtained in the previous 
calculation. This approach- was used because. the-Scannapieco model 
has two unknown. quantities, 3 2^/9 y and V^ Y . It was felt that the 
above procedure would yield an order of magnitude estimate of the 
gradient scale length required for-'ins lability and that a two-parameter 
study was not justified at this time. 

The characteristic lengths of the gradients required are on 
the order of 50-100 km, which means that the horizontal gradient of 
either n or v. must be of the same order of magnitude as the vertical 
gradient of these two quantities. This is. not an unphysical requirement, 
but is much larger than is normally seen at Arecibo. 

It should be noted that the gradients required to give growth 
times of 30 minutes are larger than those required to give growth 
times of 10 minutes because V.^ , - which was calculated from the 
Perkins model, is smaller. 

The only/ information that can be obtained from Tables 1-4 is 
that the Perkins model and the Scannapieco model are capable of 
predicting reasonable growth times under the proper set of conditions. 
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Unfortunately, there are not enough data, to-describe the full set of 
conditions that exist-during the night. A judgement about whether 
the models actually can predict F g will have to wait until simultaneous' 
measurements are made 'of the east-west neutral winds, the- 
horizontal gradients, and the drift velocities. The direction of 
propagation of the irregularities, although ignored in this thesis by- 
assuming;k is in a direction which yields maximum growth, is also 
critical to both the Perkins and Scannapieco models. 

4.2' Extension of the Perkins. Model 

The expressions for the growth rates predicted by the Perkins 
model and the Scannapieco model do not include short wavelength 
damping expected from diffusion, nor do they include contributions 
from recombination or E region coupling. The effect of recombination, 
while important in the description of the equilibrium, is a trivial 
extension- to the models and.merely adds the term -[3 to the growth 
rate. Its importance will be discussed later. The E region coupling 
issue is associated with the assumption that the integrated E region 
Pedersen conductivity is small compared to the- integrated F region 
Pedersen. conductivity. As will be shown, this is frequently not the 
case. Consequently, it becomes necessary to examine when E 
region coupling may be important, and its effect on growth. Focus . 
is on the Perkins model, although much of the analysis is also common 
to the Scannapieco model. 
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4.2.1 - Short Wavelength Diffusion 

Intuitively, one would expect that the instability would be 

damped at short wavelengths by a diffusion process.. It will be- - 

demonstrated that the short wavelength damping is contained in the 

2 2 

Hall conductivity terms, which are the terms of order v.^/co^ , - 
neglected, in the - Perkins model. The Hall conductivity is: 
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Above- 150 km, v <<co -, so Equation 4. 1 can be simplified: 
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In the F region, v - n < - < to^, so the expression can be further 
simplified: 
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To extend the Perkins model to include the Hall conductivity 
terms, a new moment must be defined: 
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It is seen by comparison with Equation 4.3 that ’the quantity 

4 1 Is the F region integrated Hall conductivity. 

The moment equations of the Perkin's model including the-’ 

• - - 2 2 

terms of order (V. /co in Equation 2.4 are: 

in ci u 
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Since 4V Sis of the order ( v ^ n / ;to c j|).> terms containing 4 J are 
negligible when compared with 2 terms , of -similar structure. This 

g ' g ^ 

means, for example, that ■ q - j — "dx. ’ su SgP s ting.'that- the;, gradients- 
involved, whether equilibrium gradients or gradients from the 
perturbation, tend, to have similar scale lengths. Equations 4. 5-and 
4. 6 can now be reduced: 
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Using 4V 2 v i n '/ a) c j^ one can obtain an es timate of the wavelength 
of the perturbation required for this term to be important: ' 


/ — * 



/-V / ’ 


65 meters 


_ O 

where the following parameters were used: = 10” , H = 50 km, 

kT = 10 erg, co . = 200 sec and t§~ E = 20 m/s. 

Cl ±5 x\.JiiQ Y 

2 

To find the contribution of the growth rate from the V . 4* term, 

1 

an expression for 4> must be obtained. It is possible to derive a 
fourth moment equation describing the time evolution of the Hall 

conductivity by multiplying the ion continuity equation by -g- (v.^/ 

2 ’ ' 
and integrating along a field line. 
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in ci ' ' 

One approach for obtaining 4 1 would be to use Equation 4. 9o 

This appears to be straightforward if one neglects terms of order 
3 3 

v.^/oo .. Unfortunately, in the short wavelength limit, terms of 
3 3 

order v. /w . can become important. One is then not able to close 
m ci r 

3 3 

the set of moment equations, since terms of' order vf /u . involve 

in ci 

higher order moments. 
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The effect of this truncation problem is only evident at 
short wavelengths. If one carried the moment equations to any order, 
and .truncated them to- that order, the growth rate would be the 
Perkins growth rate {Equation- 2. 14) for all but the short wavelengths, 
where the expansion method would break down. . 

It is clear that a method to close the moment equations must 
be adopted that does not involve arbitrary truncation, if the short 
wavelength dependence is desired. Two ways will now be presented 
that give essentially the same results. 

The first method is to express 4 1 in terms of £, in an ad-hoc 
manner. Based on- the definition of 4* , the following seems a good 
choice: 
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Using the definition of 2 and N, one obtains: 
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(4. 10) 


If it is assumed that N is constant during growth, equations 
4, 7 and 4. 8 can be linearized using: 
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It is easily shown that for perturbations of the form 

1.1c 9 -J- V ^ ^ 

e — — ip < t the term V i[; adds a damping term to the growth 

rate: 
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The additional term will dominate as the wavelength approaches, 
values on the order 20-40 meters. 

The equilibrium solution is the same as was found for the 
Perkins model,, There is no reason to suspect gradients in the Hall 


conductivity in the absence of gradients in the Pedersen conductivity. 
Thus, Equation 2.12 still describes, the equilibrium condition, 

-A somewhat "ad-hoc 1 . 1 method was adopted to close the moment 
equations, so an alternative method, .based on the use of Equation 4., 9 
is now presented. Assuming that the changes in 4 1 in time and space 
are small in the equilibrium, one obtains: 
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, and the mass of the neutrals is assumed to be 


the same as the mass the ions in the F region, one can write: 
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Assuming that the perturbed state does not deviate greatly 
from the equilibrium, 4 , q , and N Q are replaced by 4b 2, and N. 
The result: 


4 


ec 


N 


differs from Equation 4.10 by. a factor of 3, so the damping contribution 
to the growth rate will only differ by a factor of This would enable 

the perturbations of wavelengths shorter than 35 to 70 meters to be 
damped. 

It has been shown in this section that there are problems 
associated with the closure of the Perkins moment equations when the 
wavelengths of the perturbed plasma-layer becomes short. It is 
necessary to adopt a closure scheme based on intuition, and extreme 
care must be taken when the equations are truncated to any order. 

It has been, demonstrated that the short wavelength damping 
term is contained in the Hall conductivity terms that were neglected 
by. Perkins {1973). Two methods of closure were presented that gave 
essentially the same results, except for a numerical coefficient. 

Neither of these approaches is fully satisfactory, so the question of 
the correct numerical coefficient for. the short wavelength damping 
term can not be considered resolved at this time. 

4.2.2 E Region Coupling 

.The Perkins model assumes that the E region contribution to 
the field integrated Pedersen conductivity is negligible. However, 
Zinchenko (1976) found that at the Arecibo Observatory the E and F 
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region Pedersen conductivities can be of the same order magnitude. 

t 

(Figures 66 and 67). ijle also found that the E region Hall conductivity- 
can be a factor of 10 greater than the F region Pedersen conductivity, 
but this does not turn out to be of significance in the present work. 
This point will be discussed later. 

In this section, tile effect of the large E region conductivities 
will be discussed in the framework of the flux tube model of Perkins 
(1973). 

It is first noted from Figure 7 that the E and F regions can be 

treated as two separate regions. For the Arecibo data, ornwhich the 

Zinchenko (1976) calculations are based, the E region extends from 

82 km to 180 km, and the F region is considered to be heights greater 

than. 180 km. One can not, for the Perkins model, extend the limit 

of integration into the E region because in the E region v.^/co . I s 

of order unity (or larger). An expansion of c- and cr in powers of 

Jr ri 

v. /w . is not valid, 
in ci 

Because of this, it is not possible to develop a closed set of 
moment equations for the E region by the same method as that used 
by Perkins (1973). However, it is still possible to develop a model 
based on physical intuition which includes the effects of the high E 
region conductivities. This model can be explained with the help 
of Figure 68 . 

As shown, the E and F regions are two separate regions of the 

ionosphere connected by the magnetic field lines. The resistance to 

electron flow along the field lines is almost zero, so electrons can 

easily flow between the two regions. The Perkins model of F states 

s 
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Figure 66; Time va.ria.tion of integrated IT region Pedersen 
•conductivity, and E region Pedersen. and-Hall 
conductivities (from Zinchenko, 1976). 
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Figure 67: Time variation of integrated F region Pedersen 
conductivity, and E region Pedersen and Hall 
conductivities (from Zinchenko, 1976). 
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that a discontinuous electric field will be produced when the F region 

is perturbed in. the manner shown in Figure 4. 

« 

Becaus.e of .the, high .conducti-v-ity- along- -B, the field lines are 
assumed to be equipotenfcials . This will cause the same''dis continuous 
field to appear in the E region. If the E region has a higher 
transverse conductivity than the F region, the particles can flow 
across field .lines easily in. the E region, (Path II) and equalize the 
electric fields in the E region. Because the field lines are equipotentials, 
electrons will then flow up the field lines to the-F region and equalize 
the electric fields there. Thus, the E region acts as a short circuit. 

The short circuiting currents are Pedersen currents; the Hall current 
of the E region does not enter into the model because it is perpendicular 
to. the electric field. 

In reality, there is resistance between. the E and F regions, 
Zinchenko and Nisbet (1976) included the resistance between, the two 
regions and found that the E and F regions can- be considered de- coupled 
when perturbations with wavelengths less than a few kilometers are 
considered. For longer 'wavelength perturbations, the two regions are 
coupled. Dyson et al„, (1974) have observed that the irregularities 
associated with F g do not usually have wavelengths longer than 10 km. 

It is suggested that the-E region short circuiting mechanism might 
be the cause. 

The two region model shown. in Figure 68- is developed • 
mathematically by first assuming that the E region conductivities 
are not perturbed by the instability developing in the F region.- 
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Figure 68: E region coupling model 
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The assumption of plasma quasi- neutrality including the 
E region now states that: 

* (jf + ) = 0 (4, 14) 

The major currents are assumed to be the E region- Hall and Pedersen 
.currents, and the F region Pedersen, current. 
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gradients in the condufctivities 

in the E region, 

Equations 

4.14, 4, 15 and 4. 16- can be 

combined to yield the equation: 
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As was expected, the Hall current does not enter the analysis. 

The analysis in Chapter III demonstrated that for the Perkins 
instability (horizontal uniformity in S q ), the effects of recombination 
and the resultant electric field adequately model the .time evolution 
of the F region Pedersen conductivity. 
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By linearizing Equation 4. 18, and assuming. a solution of the form 
e — — + the growth rate is found to be: 
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The next step is to linearize Equation 4. 17, solve for 
and substitute the result in Equation 4,19. This yields: 
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By using several trigonometric identities, >one can rewrite the above 
expression as: 
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The angles Q and 0 are defined, in Figure 3. 

Maximum growth occurs at a = 9 / 2. 
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(4. 20) 


The short wavelength dependence has been added for completeness. 

Equation 4. 20. is the final form of the growth-rate for a model 
that includes recombination, E region- coupling, and short wavelength 
diffusion. The first term is readily- identified as the growth rate of 
the Perkins model. The additional terms all contribute to damping 
of the instability and. will be important under different conditions. 

The term describing the .wavelength dependence of the growth rate- is 
only significant at short wavelengths. Assuming that T = 1000°K 
and S^q/N^q = 0. 1, which are over estimations so. that the maximum 
damping. effect can be realized, the wavelength term can be written as;- 



This term would have to be on the order of 1/600 sec -1 to- 1/1800 sec' 
to damp the positive term from the Perkins growth rate calculated 
in- Tables 1-4. To yield this result, the -wave length would have to be 
, on- the order of 25 to 40 meters. This is normally much. less, than the 
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wavelength, of the irregularities associated with F , (Dyson et al. , 

1974), suggesting that a damping mechanism may well exist. 

The recombination term is only important when the Layer is 

low. Observations indicate that F usually does not occur when the 

s 

layer is low; damping by recombination may be a contributory factor 
under those conditions. 

For the September night, (3 = 1.47x10 ^ sec ^ at'22.47 when 

- 5 - 1 

the layer was high, 1.69 x 10 sec when onset occurred around 

-4 - 1 

21.83, and 5.4 x 10 sec at 0400 when the layerrwas low. The 

only significant addition to the growth rate for these three cases is 

at 0400. Assuming that the positive piece, which is the Perkins 

growth rate, is in the range of 1/600 sec ^ to 1/1800 sec , 

- 4 - 1 

recombination will subtract values of the order of 5.4 x 10 sec 

-4 - 1 

from this at 0400. Because 5.4 x 10 sec is around 1/1800 sec, 

recombination would obviously have an effect, especially for the slower 

growth' rate. In general, however, for onset conditions the layer is 

high enough that recombination will be a negligible effect. 

The E region coupling term is hard to analyze because of the 

degree of uncertainty in E region conductivities. The source of the 

uncertainty becomes clear if one looks at the amount of fluctuations 

in the density profiles in the E region shown in Figure 7. However, 

if one looks at integrated conductivities only, as Zinchenko (1976) 

did, the integration process smooths some of the fluctuations. 

E 

Values of S n and S are taken from Figures 64 and 65 

ir - J? 

at the onset times. These. are: 



160 


September 17: 

E 

Sp = . 1 -mhos 


S „ = .05 mhos 
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May 18: 
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S„ = 1 mho 

If 



,E 


+ 2 . 


F0 


. 17 


Equation 4. 20 is rewritten assuming that the layer- is high, so 
recombination can be neglected and the wavelength of the perturbation 
is long so the diffusion term can, be neglected: 
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The values of calculated in Tables 1-4 were such that growth 

rates of 1/600 sec ^ and 1/1800 sec ^ were obtained from the first 
piece of this expression, which will be called Y.p- 
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where y is in the range of. 1/600 sec to 1/1800 sec . To produce 
the smallest E region effect,, set 9 equal to ir/2- e, , where e is just 
large enough to yield an equilibrium eastward field. - Assuming that 
e is small, one can write: 
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Thus, for the September night at onset* the E- region coupling will 
reduce. the growth rate by at least 6*7 .per cent, while for.the May 
night,, it is reduced by. only 17 per cent. Clearly, the'E region, 
coupling can have- a significant effect on the .growth -rate., 



CHAPTER V 


SUMMARY AND CONCLUSIONS 

In this, thesis the current models of mid- latitude _F g were 
studied in detail.- The Perkins flux tube model, which has led to two 
different mechanisms for F region' instability, depending on the 
equilibrium state, seems especially promising for describing the 
nighttime F region at midlatitudes . Because particles are tied 
tightly to the field lines in the F region, the concept of a flux tube is 
a physical one. Also, the premise of the model that the nighttime 
ionosphere is supported by Pedersen currents and neutral winds 
agrees with the findings of Behnke and Kohl (1974). 

To -compare the predictions of equilibrium theory with data 
from the Arecibo Observatory, the Perkins model was then extended 
to include recombination and boundary terms. Also, the neutral wind 
was added explicitly to the model by deriving the equations in a 
stationary reference frame. 

Data from four nights of observations at the Arecibo Observatory 
were analyzed. It was found that the resultant eastward electric -field, 
defined as 

B E RESY = + ~B E oy * V NX 

is extremely well correlated with the motion of the F layer. This 
indicates that Perkins (1973) was correct in assuming that the layer 
is supported by the eastward electric field and southward neutral 
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wind (of. Equation 2. 12). Good correlation was also found between 
the field integrated conductivity and the motion of the layer., also 
supporting "the Perkins model,,. 

The flux data were immediately seen to be in poor agreement 
with the rest of the data. Above 300 km there should be no significant 
change .in the flux with altitude, because there’ are no sources or, 
losses. The data ‘showed large differences. It was found that the 
flux data did not agree at all with the data on the integrated plasma 
content. It was first thought that. perhaps the interval of integration 
was too large, thus accentuating the difference between the top flux 
and the bottom flux. Also, there were times that the S/N ratio was 
poor at the boundaries’ because the layer moved toward-the edge of 
the integration interval. 

■ A method of analysis was proposed, in which the interval of 
integration was narrowed and. followed the vertical motion of the peak. 
This guarantees a good S/N ratio at the boundaries, • Nonetheless-, 
this method of analysis did not resolve the difficulties associated 
with the flux data because the flux difference for even adjacent heights 
was too large. The net flux into or. out of a flux tube is much tod 
large and turbulent,- while the total. content of the flux tube does not 
change much. There are only two ways to explain this flux behavior: 
either there, is a large density gradient moving .into the region of 
observation, or the. flux data contains much greater errors than has. 
previously been assumed. It was concluded tha.t the latter, is the more 
likely possibility. 
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The Perkins' model was found to predict the general .time 
behavior of the plasma content and Pedersen conductivity of a flux 
tube when recombination is added. Recombination was' found to be 
a necessary addition to the model when the layer is low. 

The qua si- equilibrium of the Perkins model was compared with 
the data. It was found that the modeLagreed with the data for a large 
percentage of the time for three of the four days studied. In comparing 
the days withF g , it was found that the quasi- equilibrium assumption 
around the .time of onset was -a- good one for only one of the days. For 
the other day, the. time scale for changes in S -andJE^ were much too ■ 
rapid to assume a quasi- equilibrium. 

An attempt was made to analyze, the growth-rates predicted 
by the Perkins model and the Scannapieco model using the Arecibo 
data. However,, it was immediately found that there are not sufficient 
data. to calculate even linear growth -rates . The quantities that a're 
missing are the direction of propagation of the irregularities, the 
east-west component of the neutral wind, and the horizontal gradient 
of the field integrated Pedersen conductivity. To properly analyze the 
models of F , an experiment must be designed expressly for that 
purpose. The following quantities must be obtained simultaneously: 

(1) drift velocities; .(2) east- west neutral winds; (3) density profiles; 
(4) propagation direction of the irregularities and wavelength/if 
possible; and (5) horizontal gradients of the field integrated Pedersen 
conductivities. Only when all of these quantities are obtained - 
simultaneously will complete judgement of all the models be possible. 

It was found, however, that the dnagnitude of the neutral wind velocity 
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required to give reasonable growth times in the Perkins model, is 
plausible. Also, the Scannapieco model will predict reasonable 
growth times if the character is tic length of the east-west gradient, in 
•the Peder.sen conductivity' is on the, order of 100 km. 

The growth' rate formula of the .Perkins model was extended 
also. -It was demonstrated that a short wavelength damping 'term is 
contained in the Hall conductivity terms that were neglected by< 

Perkins .(1973). Recombination was added as a damping term that is 
only, important when the layer is low. A model was also developed 
to include the effects of the high E region conductivity. It was shown 
that the E region can short-circuit the Perkins -instability for long 
wavelengths (at least a few kilometers), if the E region conductivities 
are large enough. 

Further study of mid- latitude F is clearly required. An 

s 

experiment should be designed solely for this purpose. In the . 
meantime, there is much work to be done on ■ the data that is presently . 
available. A detailed study of more' nights should be done to determine 
the qualitative behavior of ionospheric parameters both during F 

s 

conditions and quiet conditions. In addition, data from other mid-latii- 
tude stations should be analyzed. Density contour maps, which can 
be obtained from the available Arecibo data, would help resolve the 
question of the magnitude of horizontal density gradients. Clearer 
ionograms should be obtained; with higher resolution ionograms the 
intensity of F vs. the equilibrium ionospheric parameters can be 
studied. It is also possible to obtain information. about the wavelength 
of the irregularities from good ionograms. 
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The model of E and'F region coupling proposed in this thesis 
.is clearly, an oversimplification. A better model would include the 
resistivity between the E and F regions. Also, the response of the 
E region to F region instabilities should be included, rather than 
tr eating' *the 'E region as a constant sink, as was done in the proposed 
model. 

This thesis has concentrated on equilibrium parameters and 
.linear instabilities. The turbulent structure is also important. To 
study this, a non-linear computer simulation of the Perkins instability 
-is required, similar to that performed by Scannapieco et al., (1975). 
Short pulse radar techniques (Rino and Livingston, 1974) could 
possibly supply information on the turbulent structure to compare 
with the simulations. 

Finally, the question of flux behavior-requires- further investi- 
gation. Contour plots of the plasma velocities and of the flux may 
help to illuminate the flux behavior. A careful evaluation of the method 
of data reduction of plasma velocities should be undertaken with the 
view of obtaining better error estimates. 



APPENDIX 


DESCRIPTION OF ARECIBO MODEL USED FOR 
CALCULATION OF COLLISION FREQUENCY 

A base collisionYrequency at 102 km, and a temperature 
profile for the heights 100 km to 137 km must be supplied. The base 
collision frequency was chosen as 3000 sec and the temperature 
profile is 


H (km) 

T (°K) 

99. 1 

197. 0 

102.5 

200. 0 

106. 0 

212. 0 

109. 5 

241.0 

113.0 

286. 0 

116.4 

335, 0 

119.9 

376. 0 

123.4 

423. 0 

126. 9 

464.0 

130.4 

497. 0 

137. 3 

520.0 


The neutral atmosphere densities at the base are then 
calculated, 

n, = 2.68 x 10^ * (v. )' 
b ' mb 

n b (N 2 )= .759 *n b 

n b ^2 ^ = • 17 58 * n b 

n, (0) = o 0558 * n. 
b b • 
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The exospheric temperature (T gx ) is defined as the average 
of the temperatures measured at 317.6 km and 346.,5.km. The 
temperature profile above 137 km is calculated using the following 
model: 


T(h) = T - (T - T,-- n ) e "“ 028 (h_i37 ’ 0) 
' ' ex ex 137.0 


The neutral atmosphere density above the base (102 km) is 
calculated using the expression: 


^ _ ^base x "^base _-h/H 

IN — m ^ 


KX 

where H = » and M is the mass of the neutral species whose densi ty 

is being calculated. 

The collision frequencies are calculated for various height 
ranges using the expressions: 

102-215 km: 

v. n = 3.6xl0“' 10 *n (N 2 ) + 4. 368 x 10" 10 * n <0 2 ) 

+ 2. 02 x 10 -10 * n (0) 


215-463 km: 


v._ = 5. 855 x IQ" 10 * n (N 2 ) + 6.8x10 10 * n.(0 2 ) 


in 


+ 7. 3 x 10' 10 * n (0) 
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The collision frequencies below the base are calculated using: 


v. 

in 


( V in^b 


h/H 

e 


A mass of 28. 0 is used for this expression. 

The neutral atmosphere densities are extended below the base 
using the expressions: 


DENS 


N, e 
base 


n(N 2 ) • 75 9 x N dens 

n(0 2 ) = (. 1758 H- 10 2-82 ^ x N DENS 

n(O) = -0558 (1.0 - 1 [ 5 ^ 5 -)xN dens 
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Imol, George, A Comparison Between the Currant Models of Mid-Latitude Spread F 
and Data frora the Arccibo Observatory, Tha Ionosphere Research 
Laboratory, Electrical Engineering East, University Park, Pennsylvania, 
16802, 1977. 

The current nodal* of raid-latitude T a arc studied The assumption* and 
derivation of tho Reid model, the Scannapicco model, and the Perkins raodol aro 
presented in detail. 

Incoherent-* cat ter data of tho density profile* and valocity profila* 
were obtained froo tho Arccibo Observatory In order that tho cmdols could be 
evaluated on the basis of experimental data. Initial studies indicated that the 
Perkins model was most roprescntotlvo of che data from Aroclbo, so a detailed 
comparison of the predictions of tha Perkins model and the data was made. Two of 
four nights studied are night* with F s . 

The Perkins modal is derived in a frame of reference moving with tho 
velocity of the neutral wind; the model is transformed to the rest frame to 
facilitate comparison with data, Sovcral data handling technique* arc introduced. 

In particular, an Integration interval that remains constant in length, but follows 
tha vertical motion of the peak of the F layer is used to obtain the field 
integrated quantities of the Perkins model. 

1c was found that the Perkins model describes the general time behavior 
of the field integrated Pedersen conductivity if recombination is included, and 
tha boundary flux is neglected. Tho numerical agreement 1 b usually within a 
factor of two. It was found that tho flux data Is inconsistent with the rast of 
the data. Error* on the order of 30 a/a in the Individual velocity measurements 
are believed to be present to account for the differences. 

The Perkins model Is abla to predict reasonable growth rates of the 
development of F B , howevor, a direct comparison between the growth rate* predicted 
by tho model and actual onset times is not poooiblo because of tho lack of data on 
the eastward neutral wind 

tha linear growth rate of tha Perkins model Is extended by including 
recombination, and a short wavelength damping term and coupling to the E roglon. 

It 1* found that, for wavelength perturbations several kilometers or longer, the 
E region can short circuit instabilities in the F region and reduce tho growth rata. 
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imaL, George, A Cosparluon Between the Current Models of Mid-Latitude Spread F 

and Data from the Arccibo Observatory, The Ionosphere Research 

Laboratory, ElectricaL Engineering East, University Park, Pennsylvania, 

16802, 1977 

Tha current model* of mid-latitude F s arc studied The assumptions and 
derivations of the Retd raodol, the Scannapicco model, and the Parkin* modal are 
presented in detail. 

incohcrent-scattcr data of the density profile* and velocity profiles 
were obtained from the Arccibo Obsorvatory in order that tha models could be 
evaluated on the basis of experimental data. Initial studios indicated that the 
Perkins model wa* moot repraecntative of the data from Arccibo, ao a dotailed 
comparison of the predictions of the Perkins modal and the data was made Two of 
four nights studied arc nights with F„ 

The Perkins model is derived in s frame of reference moving with the 
velocity of the neutral wind, the modal is transformed to the rest frama to 
facilitate comparison with data. Several data handling techniques arc introduced. 

In particular, an integration intorval that remains constant In length, but follows 
the vertical motion of the peak of the F layer la used to obtain the field 
integrated quantities of the Pork ins model. 

It va* found that tho Porkin* model describes tho general time behavior 
of tho field integrated Pederson conductivity if recombination is included, and 
tho boundary flux la neglected. The numerical agreement i* usuolly within a 
factor of two. It was found that tho flux data la inconsistent with the rest of 
the data. Error* on the order of 30 m la in the individual velocity measurements 
arc believed to be present to account for tho differences 

The Parkins modal is able to predict reasonable growth rates of tho 
development of F # , however, a diroct comparison between tho growth ratos predicted 
by the model and actual onset times is not possible because of tho lack of data on 
the eastward neutral wind. 

The linear growth rate of the Perkins modal is extended by Including 
recombination, and a short wavelength damping term and coupling to the E region 
It is found that, for wavelength perturbations several kilometers or longor, the 
E region can short circuit instabilities In the F roglon and reduco the growth rate. 
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Imcl, George, A Comparison Between tha Current Models of Mid-Latitude Spraad F 
and Data from tho Aroclbo Observatory, The Ionosphere Research 
Labontoty, Electrical Engineering East, University Park, Pennsylvania, 
16802, 1977 

Tho current models of raid-latitude F s aro studied. The assumptions and 
derivations of the Reid model, tho Scannapicco model, and the Perkins model arc 
presented In detail. 

Incohorcnt-s carter data of the density profiles and velocity profiles 
wore obtained from the Arccibo Observatory In order that the models could be 
evaluated on the heal* of experimental data. Initial studies Indicated chat the 
Perkins model wa* most representative of the data from Arccibo, so a dotailed 
comparison of tho predictions of Che Perkins modal and the data was made. Two of 
four nights studied are nights with F« 

Tho Perkins model la derived in a frame of reference moving with tha 
velocity of the ncutrol wind, tha model 1* transformed to the rest frame to 
facilitate comparison with data Several data handling techniques are introduced 
In particular, an Integration Interval that remains constant in length, but follows 
the vertical motion of tho peak of tha F layer la used to obtain the field 
lntcgtated quantities of the Porkin* raodol. 

'It was found that tho Porkin* model doocribes the general time behavior 
of the field integrated Pedersen conductivity if recombination is included, and 
tho boundary flux is neglected. Tho numerical agreement is usually within a 
factor of two It wo* found that tho flux data is inconsistent with the rest of 
the data. Errors on tho order of 30 m/s in the individual velocity measurements 
ore believed to be present to account for tho differences. 

Tho Perkins model is able to predict reasonable growth ratos of tho 
development of F >r however, a diroct comparison between tho growth rates predicted 
by the raodol and actual onset times is not possible because of tha lack of data on 
the eastward neutral wind. 

The linear growth rate of Che Perkin* model is extended b> including 
recombination, and a short wavelength damping terra and coupling to the E region 
It is found that, for wavelength perturbations several kilometers or longer, tho 
B region can short circuit instabilities in the F roglon and reduce the growth rate. 



1.5.3 F-Region 


-I 


Imel, CcOrga, A Comparison Between tha Current Models of Mid-Latitude Spread F 
and Data from tho Arccibo Observatory, The Ionosphere Research 
Laboratory, Electrical Engineering East, University Park, Pennsylvania, 
16802, 1977. 

Tho current models of mid-latitude P 9 aro studied The assumption!! and 
derivations of the Reid model, the Scannapicco modal, and the Perkins model are 
presented in detail 

Incoherent-scatter datn of tha density profiles and velocity profiles 
were obtained from the Arccibo Observatory in order that the models could bo 
evaluated on the basis of experimental data. Initial studies indicated chat the 
Perkins model was moot representative of the data from Aroclbo, so a detailed 
comparison of the prediction* of the Perkin* nodal and cho data wa* made. Two of 
four night* *tudlcd arc nights with F s . 

The Perkins modol ia derived in a frame of reference moving with the 
velocity of the neutral wind, the nodcl is transformed to tho rest frano to 
facilitate comparison with data. Several data handling techniques aro introduced. 
In particular, an integration intorval that remain* constant in length, but follows 
the vertical motion of the peak of the F layer is used to obtain the field 
integrated quantities of the Perkins model. 

It was found that the Perkins model describes the general time behavior 
of the field integrated Pedersen conductivity if recombination is included, and 
the boundary flux is neglected The numerical agreement is uaually within a 
factor of two. It van found that tho flux data i* inconsistent with the rest of 
tho data Errors on the order of 30 m/s in the individual velocity measurements 
arc believed to be present to account for tho differences 

The Perkins nodal is able to ptedict reasonable growth rates of tha 
development of P s ; however, a direct comparison between the growth rates predicted 
by the model and actual onset times is not possible because of the lack of data on 
the eastward neutral wind. 

The linear growth rata of tha Parkin* model is axeended by including 
recombination, and a short wavelength damping term and coupling to the E region. 

It is found that, for wavelength perturbations several kilometers or longer, tho 
E region can short circuit instabilities in tha F region and reduce the growth rate 
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I Baal, George, A Comparison Between the Curront Model* of Mid-Latitude Spread 7 

and Data from the Areclbo Observatory, The Ionosphere Research » 

Laboratory, Electrical Engineering Ease, University Park, Pennsylvania, 
16802, 1977. 

The current oodcln of aid-la tituda F s are studied. The assumptions and 
derivations of the Reid model, the Scsnnapiaco model, and the Perkins model arc 
presented in detail. 

Incohcrent-Bcatter data of the density profiles and velocity profiles 
vc re obtained from the Areclbo Observatory In order that the models could be 
evaluated on the basis of experimental data, Initinl studies indicated that the 
Parkins modal vas most representative of the data from Areclbo, so a detailed 
comparison of the predictions of the Perkins model and the data was made. Two of 
four nights studied aro nights with P_. 

The Perkins modeL Is derived in a frame of reference moving with tho 
velocity of tho neutral wind, the model is transformed to the rest frame to 
facilitate comparison with data Several data handling techniques arc introduced 
In particular, an integration interval that remains constant in length, but follows 
the vertical motion of the peak of the F layer Is used to obtain tha field 
Integrated quantities of tha Perkins model. 

It vas found chat the Porklns model describes the general time behavior 
of the field integrated Pedersen conductivity if recombination is included, and 
the boundary flux la neglected. The nunorlcal agreement is usually within a 
factor of two It van found chat tho flux data is Inconsistent vlth tho rest of 
the data. Errors on Che order of JO a/s in tha individual velocity measurements 
aro believed to bo present to account for the differences. 

The Parkins modal is able to predict reasonable growth rates of tha 
development of F s ; however, a direct comparison between the growth rates predicted 
by the model and actual onset times is not possible because of the lack of data on 
the eastward neutral wind. 

The linear growth rate of tho Parkins modal is extended by including 
recombination, and a short wavelength damping term and coupling to the E region 
It i* found that, for wavelength perturbations several kilometers or longer, the 
E region can short circuit Instabilities in the F region and reduce the growth rate. 
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Imol, George, A Comparison Between tho Currant Models of Mid-Latitudo Spread P 
and Data from tho Areclbo Observatory, The Ionosphere Research 
Laboratory, Electrical Engineering East, University Park, Pennsylvania, 
16802, 1977 

The current models of mid-latitude P s are studied. The assumptions and 

derivations of the Roid model, tha Sconnaplaco model, and tho Porklns model aro 
presented In detail. 

Incohc rent-oca ttcr data of tho density profiles and velocity profiles 
were obtained from tho Areclbo Observatory In order that tho models could be 
evaluated on the basis of experimental data. Initial studios indicated that the 
Perkins model was most representative of the data from Arocibo, so a detailed 
comparison of the predictions of tho Perkin* nodal and the data vaa nado. Two of 
four nights studied arc nights with ? 0 . 

The Perkin* nodal is derived in a frame of reference moving with the 
velocity of the neutral wind, the model is transformed to the rest frame to 
facilitate comparison with data Several data handling techniques arc introduced* 

In particular, an Integration interval that remains constant in length, but follows 
tho vertical motion of the peak of the F layer is used to obtain the field 
integrated quantities of the Perkins modol. 

It was found that the Perkins modal describes tha general time behavior 
of the field integrated Pedersen conductivity if recombination is included, and 
the boundary flux Is neglected. The numerical agreement is usually within a 
factor of two. It wae found that the flux data is inconsistent with the rest of 
the data Errors on the order of 30 m/s in the individual velocity measurements 
are believed to be present to account for the differences 

Tha Perkins model is able to predict reasonable growth rates of the 
development of P fi> however, a direct comparison between the growth rates predicted 
by the model and actual onset clmos is not possible because of the lack of data on 
the eastward neutral wind. 

Tha linear growth rate of the Perkins model is extended by including 
recombination, and a short wavelength damping term and coupling to the E region. 

It is found that, for wavelength perturbations several kilometers or longer, the 
E region can short circuit instabilities in tho F roglon and reduce the growth rata. 
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Xmel, George, A Comparison Between tho Curront Models of Mid-Latitude Spread F 
and Data from the Areclbo Observatory, The Ionosphere Research 
Laboratory, Electrical Engineering East, University Park, Pennsylvania, 
16802, 1977. 

Tho current models of mid-latitude F c aro studied. The assumptions and 
derivations of tho Reid model, the Scannapieco model, and the Perkins model are 
presented in detail. 

Incohorent-scattor data of the density profiles and velocity profiles 
were obtained from the Areclbo Observatory in order that tho models could be 
evaluated on the basis of experimental data, IalClal studies Indicated that the 
Perkins model was most representative of the data from Areclbo, so a detailed 
comparison of Che predictions of the Porklns model and the data was made Two of 
four nights studied are nights with F s . 

The Perkins model lc derived in a frame of reference moving with tho 
velocity of the neutral wind, tho model is Cron* formed to tho rest frame to 
facilitate comparison with data Several data handling technique* arc introduced. 

In particular, an integration interval that remains constant in length, but follows 
the vertical motion of the peak of the 7^ layer is usod to obtain the field 
Integrated quantities of the Parkins modal. 

It was found chat the Perkins model describes tho general timo behavior 
of the field integrated Pedersen conductivity if recombination is included, and 
the boundary flux is neglected. The numerical agreement is usually within a 
factor of two. It was found chat the flux data is inconsistent with the rest of 
the data. Errors on the order of 30 ra/s in the individual velocity measurement* 
are beliovod to be present to account for the differences. 

Tha Perkins model is abla to predict reasonable growth rote* of the 
development of F 0 , however, a direct comparison between the growth rates predicted 
by the model and actual onset times 1* not possible because of tho lack of data on 
the eastward neutral wind. 

The linear growth rata of tha Perkins model is extended by including 
recombination, and a short wavelength damping tarn and coupling to Che E region. 

It is found that, for wavelength perturbations several kilometers or longer, the 
E region can short circuit Instabilities in the F region and reduce tho growth rato. 
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I axil, George, A Comparison Between the Current Models of Mid-Latitudo Spread F 
and Data from the Areclbo Observatory, The Ionosphere Research 
Laboratory, Electrical Engineering East, University Park, Pennsylvania, 
16502 , 1977. 

the current models of mid-latitude P 8 are studied. The assumptions and 
derivations of the Reid model* the Scannapieco model, and the Perkins model are 
presented In detail. 

Incohcrcnt-acattcr data of tho density profiles and velocity profiles 
were obtained from the Areclbo Observatory in order thae the module could bo 
evaluated on tho baoia of experimental data. Initial studios indicated that tho 
Parkin# model was moat representative of tho data from Arocibo, so a detailed 
comparison of tho predictions of the Perkins model and the data was made. Two of 
four nights studied aro nights with f s 

Tho Perkins model is derived in a frame of reference moving with tho 
velocity of the neutral wind; Che model is transformed to tho rest frame to 
facilitate comparison with data. Several data handling techniques are introduced. 

In particular, an integration interval thee remains constant in length, but follows 
tlio vortical mocion'of Che peak of tho F layer la usod to obtain the field 
Integrated quantities of tho Perkin* model, 

.It vas found chat the Porklns Dodel describes Che general time behavior 
of tho field integrated Pedersen conductivity if recombination is included, and 
tho boundary flux is neglected. Tho numerical agreement is usually within a 
factor of two It was found that the flux data Is inconsistent with the rest of 
the data. Errors on Che order of 30 a/s in the individual velocity measurements 
are beliovod to be present to account for tha differences. 

Tho Porklns model Is able to predict reasonable growth ratos of the 
development of F ft , however, a direct comparison between the growth rates predicted 
by the modol and actual onset times is not possible because of the lack of data on 
the eastward neutral wind. 

The llnoar growth rate of tho Perkins model It extended by including 
recombination, and a short wavolongth damping torn and coupling to the E roglon. 

Ic is found chat, for wavelength perturbations several kilometers or longer, the 
E region can short circuit instabilities in the F region and reduce the growth rato. 
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